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1.  INTRODUCTION 


1.1  Project  Objectives 

The  objective  of  this  two-year  study  is  to  simulate  detection  and  location  capability  of 
seismic  networks  including  regional  stations  and  arrays  in  and  around  the  Soviet 
Union.  Three  specific  research  tasks  are: 

(1)  Enhance  and  validate  the  extended  version  of  the  Seismic  Network  Assessment 
Program  for  Detection  [ SNAP/D ,  Ciervo  et  ai,  1985]  called  SNAP/DX  [Bratt  ei 
at.,  1987a]  to  accurately  represent  the  treaty  monitoring  capability  of  seismic  net¬ 
works  including  regional  stations  and  arrays. 

(2)  Normalize  SNAP/DX  to  the  observed  performance  of  existing  stations  and 
expected  conditions  in  and  around  the  Soviet  Union. 

(3)  Apply  the  normalized  simulation  methods  to  assess  the  treaty  monitoring  capabil¬ 
ity  of  existing  and  proposed  seismic  networks. 

1.2  Current  Status 

The  effort  during  the  first  year  of  this  project  was  divided  among  SNAP/DX  enhance¬ 
ment  and  SNAP/DX  normalization.  "Tie  first  semi-annual  report  described  the  normali¬ 
zation  results  obtained  during  the  first  six  months  of  the  project  [Sereno,  1989].  This 
included  two  studies:  (1)  a  theoretical  study  of  the  sensitivity  of  Pn  geometric  spread¬ 
ing  to  the  velocity  gradient  in  the  upper  mantle,  and  (2)  an  empirical  study  of  the 
attenuation  of  Pn  and  Lg  phases  recorded  in  eastern  Kazakhstan.  We  have  continued 
to  work  on  the  normalization  since  this  first  report,  but  the  more  recent  work  has  been 
primarily  in  assembling  a  data  base  of  regional  events  recorded  by  the  NORESS  and 
ARCESS  arrays  in  Norway.  We  have  processed  data  from  nearly  100  events  recorded 
by  the  two  arrays.  This  processing  includes  interactive  picking  of  regional  phases,  cal¬ 
culation  of  Fourier  signal  and  noise  spectra,  and  the  determination  of  phase  velocity, 
azimuth,  amplitude,  and  dominant  frequency.  These  signal  spectra  will  be  inverted  for 
source  and  attenuation  parameters  near  the  beginning  of  the  second  year  of  the  con¬ 
tract,  and  these  results  (along  with  the  rest  of  the  normalization  results)  will  be 
included  in  the  third  semi-annual  report  for  this  project 

This  annual  report  describes  our  enhancements  to  the  SNAP/DX  computer  program. 
The  most  important  of  these  enhancements  is  the  introduction  of  frequency  dependence 
into  the  estimates  for  the  source,  station  noise  (for  primary  and  secondary  phases), 
attenuation,  and  array  gain.  SNAP/D  and  SNAP/DX  calculate  detection  and  location 
capability  at  fixed  frequency.  This  is  adequate  to  represent  the  capability  of  telese- 
ismic  networks  because  the  dominant  signal  frequency  is  only  weakly  dependent  on 
distance.  However,  at  regional  distances  the  frequency  of  the  maximum  signal-to- 
noise  ratio  (SNR)  may  be  strongly  distance-dependent  Therefore,  the  detection  capa- 
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bilities  at  individual  stations  in  the  network  must  be  evaluated  at  a  frequency  that 
depends  on  epicentral  distance. 

1,3  Outline  of  the  Report 

This  report  is  divided  into  five  sections  including  this  introduction.  Section  2  is  an 
overview  of  our  new  netwoik  simulation  computer  program,  NetSim .  We  start  with  a 
brief  introduction  that  describes  the  relationship  between  NetSim  and  other  computer 
programs  that  have  previously  been  used  to  assess  the  capability  of  seismic  networks 
to  monitor  underground  nuclear  explosion  testing.  Next,  we  give  functional  descrip¬ 
tions  of  the  NetSim  detection  and  location  modules.  Section  3  summarizes  the  NetSim 
input  and  output  data  and  graphics  (Appendix  A  gives  a  detailed  listing  and  description 
of  each  input  and  output  parameter).  Section  4  gives  sample  NetSim  runs  for 
frequency-dependent  and  frequency-independent  input.  Section  5  summarizes  the  main 
conclusions. 


2.  NETSIM  OVERVIEW 


2.1  Introduction 

NetSim  is  a  computer  program  to  simulate  the  detection  and  location  capability  of 
seismic  networks  including  regional  stations  and  arrays.  The  main  advantage  of  Net¬ 
Sim  compared  to  similar  computer  programs  is  the  incorporation  of  frequency  depen¬ 
dence  in  the  estimates  for  the  source,  station  noise,  attenuation,  and  array  gain.  Com¬ 
puter  programs  that  simulate  network  performance  at  fixed  frequency  include  NET- 
WORTH  [Wirth,  1977],  SNAP/D  [Ciervo  et  al.,  1985],  and  an  extended  version  of 
SNAP/D  called  SNAP/DX  [Bratt  et  al.,  1987a].  SNAP/D  allows  specification  of  roul- 
tiwave  detection  criteria;  a  feature  that  is  not  available  in  its  predecessor,  NETWORTH. 
SNAP/DX  uses  the  same  detection  module  as  SNAP/D,  but  it  has  a  new  location 
module.  This  new  module  is  based  on  the  TTAZLOC  location  program  of  Bratt  and 
Bache  [1988],  and  uses  estimates  of  both  arrival  time  and  azimuth  standard  deviation 
to  estimate  the  location  uncertainty.  SNAP/DX  offers  two  approaches  for  approximat¬ 
ing  the  effect  of  undetected  phases  on  the  location  uncertainties;  the  probability- 
weighted  approach  used  in  SNAP/D,  and  a  Monte  Carlo  approach  that  mimics  the  way 
location  uncertainties  vary  in  practice.  NetSim  uses  the  enhanced  SNAP/DX  location 
module  combined  with  an  enhanced  detection  module  that  incorporates  frequency 
dependence  in  the  propagation  and  noise  characteristics.  It  also  includes  a  new 
representation  of  noise  levels  for  secondary  phases  that  is  based  on  the  coda  of  previ¬ 
ous  arrivals. 

2.2  Detection  Module 

The  detection  module  calculates  the  signal  and  noise  amplitudes  for  each  wave,  station, 
and  epicenter.  These  amplitudes  are  used  to  determine  the  probability  of  detection  at 
individual  stations.  The  individual  station  probabilities  are  combined  to  determine  the 
probability  of  detection  for  the  network.  liie  NetSim  detection  algorithm  is  essentially 
the  same  as  that  of  SNAP/D  and  SNAP/DX.  However,  there  are  two  main  exceptions. 
First,  the  signal  and  noise  amplitudes  at  individual  stations  are  calculated  as  a  function 
of  frequency.  Second,  the  noise  for  secondary  phases  depends  on  the  amplitude  spec¬ 
trum  of  previous  arrivals  and  the  coda  decay  rate.  Therefore,  the  noise  for  secondary 
phases  depends  on  event  magnitude  and  epicentral  distance. 

The  NetSim  detection  module  has  two  options.  It  can  calculate  the  detection  threshold 
of  the  network  at  a  fixed  confidence  level  (these  thresholds  are  determined  by  varying 
the  event  size  until  the  detection  probability  of  the  network  equals  the  desired 
confidence  level),  or  it  can  calculate  the  probability  that  the  network  will  detect  an 
event  of  fixed  size.  In  either  case,  the  calculations  are  performed  for  each  epicenter  in 
a  grid  specified  by  the  user.  Hie  results  are  displayed  on  a  map,  either  as  contours  of 
event  size  (for  fixed  confidence  level),  or  as  contours  of  detection  probability  (for  fixed 
event  size). 


3 


The  first  step  in  the  detection  module  is  to  calculate  the  probability  of  detection  at 
individual  stations.  Assuming  that  both  signal  and  noise  are  log  normally  distributed, 
the  probability  that  wave  k  will  be  detected  at  station  i  from  an  event  at  epicenter  j  is: 


where. 


Pijt  =  *,4> 


log  %  -  log  %  -  log  SNR  * 

(<4* +<%*)“ 


(2.2.1) 


(2.2.2) 


and, 


Pi 

Syk 

Nijk 

SNR* 


aSijk 


°Nijk 


reliability  of  station  i  (the  fraction  of  time  that 
station  i  is  operational). 

signal  amplitude  for  wave  k  at  station  i  from 
epicenter  j. 

noise  amplitude  for  wave  k  at  station  i  from 
epicenter  j. 

signal-to-noise  ratio  necessary  for  detection  of 
wave  k  at  station  i. 

standard  deviation  of  log  signal  amplitude  for 
wave  k  at  station  i  from  epicenter  j. 
standard  deviation  of  log  noise  amplitude  for 
wave  k  at  station  i  from  epicenter  j. 


The  signal  and  noise  amplitudes  at  individual  stations  are  calculated  at  frequencies 
selected  by  the  user.  The  individual  station  probabilities  in  (2.2.1)  are  evaluated  at  the 
frequency  that  maximizes  the  signai-to-noise  ratio,  and  this  frequency  is  generally  a 
function  of  distance.  The  individual  station  probabilities  are  combined  to  determine 
the  probability  of  detection  for  the  network.  This  calculation  is  the  same  as  the 
SNAP/D  implementation  which  has  already  been  described  in  detail  by  Ciervo  et  al. 
[1985].  The  calculation  of  the  signal  and  noise  amplitudes  and  standard  deviations  are 
described  in  the  following  two  sections. 


22.1  Signal 

The  signal  amplitudes  at  individual  stations  Sp  are  calculated  using  one  of  two  equa¬ 
tions  depending  on  whether  the  jth  epicenter  is  at  a  regional  or  teleseismic  distance 
from  the  ith  station  (the  crossover  distance  is  selected  by  the  user).  For  regional  dis¬ 
tances: 


log  Sijtf)  -  log  C-,jk  +  5j&)  +  fl*(Ajp  f)  +  SRuff)  +  stacor *  (2.2.3) 

where  Cyk  is  a  constant  that  depends  on  the  physical  properties  of  the  source  and 
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receiver  media,  sjff)  is  the  log  source  spectrum  for  the  kth  wave,  is  the 

attenuation  of  the  kth  wave  as  a  function  of  epicentral  distance  Ay  and  frequency/,  and 
SRuff)  is  a  combination  of  the  local  site  response  and  array  gain  as  a  function  of  fre¬ 
quency  for  the  kth  wave  observed  at  the  ith  station.  Each  of  these  terms  are  described 
in  more  detail  below.  The  last  term  in  (2.2.3),  stacor **,  is  an  amplitude  station  correc¬ 
tion  for  the  kth  wave  at  the  ith  station. 


The  constant  Cijk  is  expressed  as  [e.g.,  Aki  and  Richards,  1980]: 


Cijk 


109 

4k  (pj  pi  vjk  v^)1 72 


(2.2.4) 


where  p,  and  p,-  are  the  densities  of  the  source  and  receiver  media,  respectively.  Simi¬ 
larly,  vjk  and  Vi*  are  the  seismic  velocities  of  the  kth  wave  at  source  and  receiver, 
respectively.  The  velocity  and  density  at  the  jth  epicenter  are  determined  from  a 
source  medium  geographical  grid,  and  the  velocity  and  density  of  the  ith  station  are 
determined  from  a  list  of  station  parameters.  The  factor  of  109  in  (2.2.4)  converts 
amplitudes  from  meters  to  nanometers. 


The  source  spectrum  is  represented  as  a  product  of  a  wave-independent  moment  spec¬ 
trum  M0jff)  and  a  wave-dependent  excitation  factor,  k*.  Specifically,  the  log  source 
spectrum  is  expressed  as: 


sjiff)  =  log  k*  +  log  M0jff)  (2.2.5) 

Both  k*  and  M0j  depend  on  the  source  medium  (for  example,  different  source  spectra 
are  used  for  explosions  in  salt,  granite,  or  tuff).  The  factor  k*  represents  the  relative 
source  excitation  of  different  wave  types  (e.g.,  typical  values  for  explosions  are  1.0  for 
P  phases  and  0.3  for  S  phases). 

The  Bi^Ay,  f)  are  frequency-dependent  amplitude  versus  distance  curves.  Separate 
curves  are  tabulated  for  each  wave  and  for  each  propagation  medium  (e.g.,  stable  or 
tectonic).  A  propagation  grid  identifies  the  geographical  locations  of  these  different 
propagation  media.  The  effective  attenuation  for  each  wave  used  in  (2.2.3)  is  a  linear 
combination  of  the  separate  B^Ay,f)  for  each  propagation  medium  between  the  jth  epi¬ 
center  and  the  ith  station.  The  path  weights  applied  to  the  separate  attenuation  curves 
are  calculated  from  the  relative  portion  of  the  great  circle  path  in  each  propagation 
medium  [Ciervo  et  al.,  1985]. 

The  local  site  response  spectrum  SRuff)  includes  effects  such  as  free  surface 
amplification,  instrument  response,  and  array  gain.  Since  the  signal-to-noise  ratio  is 
used  to  determine  detection  probabilities  at  individual  stations,  the  actual  instrument 
response  spectrum  is  not  needed.  Instead,  the  instrument  response  portion  of  SRuff)  is 
set  to  zero  for  frequencies  included  in  the  instrument  passband,  and  it  is  set  to  a  large 
negative  number  for  frequencies  outside  the  instrument  passband.  This  is  useful  for 
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representing  the  capability  of  networks  that  include  both  broad  band  and  narrow  band 
stations.  It  is  also  important  to  represent  the  capability  of  networks  that  include  both 
single  stations  and  regional  arrays.  For  single  stations,  the  array  gain  portion  of  SRytf) 
is  set  to  zero  for  all  frequencies.  However,  since  beamforming  can  be  very  effective 
in  suppressing  uncorrelated  noise  [e.g.,  Kveama ,  1989],  the  array  gain  portion  of  SRJJ) 
will  generally  be  greater  than  zero  foi  array  stations.  For  example,  if  beamforming 
results  in  V) V  noise  suppression,  and  if  the  signal  is  perfectly  coherent  across  the  array, 
then  SRuff)  =  1/2  log  N.  Since  signal  coherence  (and  the  number  of  array  elements 
used  in  beamforming)  generally  decrease  with  increasing  frequency,  SRutf)  is  expected 
to  decrease  with  frequency  for  array  stations. 

The  signal  amplitude  at  teleseismic  distances  has  a  slightly  different  representation 
than  the  signal  amplitude  at  regional  distances  (2.2.3).  For  teleseismic  distances  the 
signal  amplitude  is  expressed  as: 

log  =  log  Cijk  +  Sjtf)  +  Bk  ^(A ij,f)  +  SRjjff)  +  stacorn  +  Scorjk  +  RcorA  (2.2.6) 

There  are  two  main  differences  between  this  representation  and  the  corresponding  rela¬ 
tion  for  regional  distances.  First,  the  path-weighted  amplitude  versus  distance  curve  in 
(2.2.3)  is  replaced  with  a  reference  amplitude  versus  distance  curve  Bk^EP(Aij,f)  that  is 
used  for  all  source-receiver  paths.  This  assumes  that  strong  lateral  variations  in  the 
earth  are  confined  to  shallow  structure,  so  that  teleseismic  waves  (which  propagate  at 
greater  depth)  travel  primarily  in  a  laterally  homogeneous  medium.  The  second 
difference  is  that  the  teleseismic  representation  includes  source  and  receiver  amplitude 
correction  factors  ( Scorjk  and  Rcor £  that  depend  on  the  propagation  medium  of  the  jth 
epicenter  and  ith  station.  These  factors  are  used  to  model  the  effects  of  near-source 
and  near-receiver  geology  that  can,  for  example,  introduce  a  magnitude  bias  between 
two  receiver  or  source  sites. 

The  variance  of  the  log  signal  amplitude  for  each  wave  is  the  sum  of  the  variance  of 
the  amplitude  versus  distance  relation  and  die  station-specific  log  amplitude  vari¬ 
ance  of*.  That  is, 

Gsijk  -  GBA&jj)  +  °i *  (2.2.7) 

For  regional  distances,  the  variance  of  the  amplitude-distance  relation  is  a  linear  com¬ 
bination  of  the  variances  for  the  individual  propagation  media.  The  weights  are  deter¬ 
mined  from  the  portion  of  the  great  circle  path  between  source  and  receiver  in  each 
propagation  medium. 


222 Noise 

The  noise  spectrum  for  each  wave  and  station  is  expressed  as  a  sum  of  the  ambient 
noise  spectrum  and  (for  secondary  phases)  a  signal-generated  noise  spectrum  which 
includes  the  coda  of  the  previous  arrival.  This  signal-generated  noise  spectrum 
depends  both  on  event  size  and  epicentral  distance.  The  noise  power  spectral  density 
for  each  wave,  station,  and  epicenter  is  expressed  as: 


PSDijtf)  =  PSDaff)  +  PSDsijtff)  (2.2.8) 

where  PSDatf)  is  the  power  spectral  density  of  the  ambient  noise  at  the  ith  station,  and 
PSDsijk(f)  is  the  power  spectral  density  of  the  signal-generated  noise  for  the  kth  wave  at 
the  ith  station  from  the  jth  epicenter. 

The  frequency  content  of  the  signal-generated  noise  for  secondary  phases  is  expected 
to  be  much  different  from  the  frequency  content  of  the  ambient  noise,  and  similar  to 
the  frequency  content  of  the  previous  arrival.  Also,  the  amplitude  of  the  noise  for 
secondary  phases  depends  on  the  amplitude  of  the  previous  arrival  and  on  the  coda 
decay  rate.  The  simplest  assumption  consistent  with  these  dependencies  is  that  the 
signal-generated  noise  spectrum  for  secondary  phases  is  equal  to  the  signal  spectrum  of 
the  previous  arrival  multiplied  by  a  scaling  factor  that  depends  on  distance  (or, 
equivalently,  on  the  time  separation  between  the  secondary  phase  and  the  previous 
arrival).  This  is  the  assumption  used  by  NetSim  to  determine  noise  levels  for  secon¬ 
dary  phases.  The  power  spectral  density  of  the  signal-generated  noise  is  calculated 
from  the  amplitude  spectrum  of  the  previous  arrival  using: 

PSDsijiff)  =  L  (2.2.9) 

where  $  is  the  coda  decay  rate,  and  T{k-\)  is  the  time  window  length  used  to  calcu¬ 
late  the  signal  spectrum  of  the  previous  arrival.  This  parameterization  assumes  that  the 
coda  decay  rate  does  not  depend  on  frequency.  A  more  complicated  model  that 
includes  frequency  dependence  in  yk  could  also  be  adopted  if  it  is  found  that  this  is 
important  for  simulating  the  capability  to  detect  regional  secondary  phases. 

It  is  straight  forward  to  calculate  the  noise  amplitude  spectrum  from  the  noise  power 
spectral  density.  The  relationship  is: 


Nijtf)  =  [TkPSDijMm  (2.2.10) 

where  PSDij/tf)  is  calculated  using  (2.2.8)  and  (2.2.9),  and  Tk  is  the  time  window  length 
used  to  calculate  the  signal  spectrum  for  the  kth  wave.  The  noise  spectrum  in  (2.2.10) 
is  used  in  (2.2.1)  to  calculate  the  probability  of  detection  at  individual  stations. 

The  standard  deviation  of  the  log  noise  amplitude  is  a  function  of  the  standard  devia¬ 
tion  of  the  log  ambient  noise  at  each  station,  and  (for  secondary  phases)  the  standard 
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deviation  of  the  log  signal  amplitude  of  the  previous  wave.  Assuming  that  signal  and 
noise  are  log  normally  distributed,  the  standard  deviation  of  the  log  noise  for  the  kth 
wave  at  the  ith  station  from  an  event  at  the  jth  epicenter  is: 


aNIJk 


\pSDaff?  oL  +  PSDsg/ff)2  o?/,-(*-i)]1/2 
PSDatf)  +  PSDSijtf) 


(2.2.11) 


where  oNai  is  the  standard  deviation  of  the  log  ambient  noise  amplitude  at  the  ith  sta¬ 
tion,  and  Osijtjt- 1)  is  the  standard  deviation  of  the  log  signal  amplitude  of  the  previous 
wave.  Noie  that  as  the  signal-generated  noise  approaches  zero,  the  standard  deviation 
of  the  log  noise  for  the  kth  wave  approaches  the  standard  deviation  of  the  log  ambient 
noise. 


2.3  Location  Module 

NetSim  estimates  location  capability  by  computing  the  dimensions  of  the  hypocentral 
confidence  ellipsoid  for  each  epicenter  in  the  grid.  The  location  uncertainty  can  be 
calculated  at  fixed  event  size,  or  at  the  detection  threshold  of  the  network.  In  the  latter 
case,  the  detection  module  is  used  to  calculate  the  detection  threshold  for  each  epi¬ 
center  before  activating  the  location  module.  This  section  gives  a  brief  description  of 
the  calculations  performed  by  the  NetSim  location  module.  More  details  are  given  by 
Bratt  et  al.  [1987a]  and  Bran  and  Bache  [1988].  Section  2.3.1  describes  the 
probability-weighted  location  uncertainties  (the  method  used  by  SNAP/D)  and  Section 
2.3.2  describes  the  Monte  Carlo  approach  for  estimating  location  uncertainties. 

The  location  confidence  bounds  depend  on  the  distribution  of  the  detecting  stations,  the 
presumed  earth  structure,  ernd  a  priori  estimates  of  the  data  variance.  Effects  of  net¬ 
work  geometry  and  earth  structure  are  reflected  in  the  matrix  of  partial  derivatives  of 
the  data  with  respect  to  the  hypocentral  coordinates  and  origin  time  [Bratt  et  al., 
1987a].  Traditionally,  data  have  been  phase  arrival  times.  However,  azimuth  data 
from  arrays  or  three-component  stations  can  also  provide  significant  constraint  on  loca¬ 
tion  solutions.  In  fact,  under  conditions  where  few  data  are  available  they  permit 
event  locations  not  possible  using  arrival-time  data  alone. 

The  NetSim  location  algorithm  computes  location  uncertainties  using  both  arrival  time 
and  azimuth  partial  derivatives  for  any  seismic  phase  [Bratt  et  al.,  1987a].  If  the  data 
variances  are  assumed  to  be  known  a  priori  (an  assumption  required  for  network  simu¬ 
lation),  then  the  points  xe  on  the  p  percent  confidence  ellipsoid  for  the  solution  x  are 
obtained  from 


(*«  -  X?  V?  (xt -x)  =  Xp  m  (2.3.1) 

where  xj  is  the  p  percent  chi-squared  statistic  with  M  degrees  of  freedom  (M  is  the 
number  of  hypocentral  parameters)  and  Vx  is  an  estimate  for  the  parameter  covariance 
matrix  [e.g.,  Jordan  and.  Sverdrup,  1981;  Bratt  and  Bache,  1988].  The  arrival  time 
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and  azimuth  variances  for  each  wave  are  a  sum  of  the  station-specific  variances  and 
the  variances  that  depend  on  the  propagation  media  between  the  source  and  receiver. 
The  path  weights  are  determined  from  the  relative  fraction  of  the  great  circle  path  in 
each  medium. 

In  the  simulations  the  individual  stations  have  some  probability  (£100%)  of  detecting 
an  event,  but  in  practice  a  phase  is  either  detected  and  used  in  the  location  solution  or 
it  is  not  detected  and  cannot  be  used  to  constrain  location.  NetSim  offers  two  ways  in 
which  the  probability  of  detecting  a  given  phase  can  be  used  to  estimate  location 
uncertainties.  One  of  these  ways  is  by  weighting  the  standard  deviation  of  each  datum 
by  the  square  root  of  the  probability  of  detection  [Ciervo  et  at.,  1985].  Another  way  is 
to  compute  a  large  number  of  confidence  ellipses  at  each  epicenter,  using  subsets  of  all 
detections  selected  in  a  random  or  Monte  Carlo  fashion  based  on  the  probability  of 
detection  [Bratt  et  al .,  1987a].  These  two  methods  are  described  in  more  detail  in  the 
following  subsections. 

23.1  Probability  Weighting 

In  the  probability-weighted  approach,  the  effect  of  undetected  phases  is  approximated 
by  assuming  that  every  phase  is  detected  and  produces  data  whose  standard  error  is 
increased  by  a  factor  that  depends  on  the  probability  of  detection  [Ciervo  et  al.,  1985]. 
The  standard  deviation  of  each  datum  a  is  divided  by  the  square  root  of  the  probability 
that  the  station  in  question  will  detect  that  phase.  The  effective  standard  deviations 
(o/a/p^  are  increased  for  phases  with  low  probability  of  detection,  and  therefore  these 
phases  do  not  contribute  much  to  constraining  the  event  hypocenter. 

The  probability-weighted  approach  for  computing  average  confidence  bounds  is  most 
accurate  when  the  parameters  of  the  error  ellipse  are  not  very  sensitive  to  the  exclusion 
of  an  individual  datum  [Ciervo  et  al.,  1985].  However,  when  few  data  are  available  to 
locate  an  event  (e.g.,  for  most  events  near  the  detection  threshold),  each  datum  is  very 
important  Thus,  the  accuracy  of  the  probability-weighted  location  uncertainties 
degrades  with  decreasing  event  size.  Bratt  et  al.  [1987a]  show  that  the  probability- 
weighted  approach  increasingly  underestimates  the  Monte  Carlo  uncertainties  as  the 
event  size  decreases,  which  suggests  that  this  method  does  not  adequately  degrade  the 
contribution  of  data  with  low  detection  probabilities. 

232  Monte  Carlo 

The  Monte  Carlo  approach  does  not  assume  tha*  every  phase  is  detected.  Instead,  a 
large  numt'er  of  confidence  ellipses  are  calculated  at  each  epicenter,  and  each  realiza¬ 
tion  includes  data  from  a  subset  of  all  phases.  The  subset  is  determined  by  comparing 
the  probability  of  detection  at  each  station  to  a  uniformly  distributed  random  number 
between  0  and  1.  If  the  probability  of  detection  is  greater  than  the  random  number, 
then  the  phase  is  included  in  the  location  confidence  calculation.  Otherwise,  the  phase 
is  not  detected  and  cannot  be  used  to  constrain  location.  From  these  realizations. 


statistics  such  as  the  median  and  90th  percentile  confidence  bounds  can  be  determined. 
7s  approach  mimics  the  process  by  which  average  confidence  ellipses  would  be 
obtained  from  many  located  events  in  a  region,  but  it  is  more  computationally  inten¬ 
sive  than  the  probability-weighted  approach. 


3.  NETSIM  INPUT/OUTPUT 


This  section  gives  a  brief  description  of  the  NetSim  input  and  output  data,  and  the 
graphics  capabilities.  Appendix  A  lists  the  contents  of  each  input  and  output  data  file 
and  gives  a  more  detailed  description  of  each  of  die  parameters. 

Figure  3.1  is  a  schematic  top-level  data  flow  diagram  for  NetSim.  The  input  data  are 
organized  into  five  categories: 

(1)  Control  Data 

(2)  Source  Data 

(3)  Propagation  Data 

(4)  Site/Station  Data 

(5)  Noise  Data 

The  data  in  each  of  these  categories  are  described  in  Section  3.1.  In  general,  each 
category  includes  several  input  files  (indicated  by  the  smallest  boxes  in  Figure  3.1). 
There  are  a  total  of  13  different  types  of  input  files.  Some  of  these  files  are  indices 
(or  pointers)  and  others  are  one-dimensional  or  two-dimensional  data  tables.  Section  4 
of  this  report  gives  examples  of  the  input  files  for  a  few  sample  runs. 

The  output  of  the  NetSim  detection  module  is  either  detection  thresholds  for  a  fixed 
confidence  level,  or  detection  probability  for  a  fixed  event  size.  Hie  location  module 
outputs  the  parameters  of  the  hypocentral  confidence  ellipsoid  for  a  fixed  confidence 
level.  The  NetSim  output  data  are  described  in  Section  3.2.  A  graphics  package  called 
SnapCon  is  used  to  contour  the  NetSim  output  on  regional  or  world  maps  [Bratt  et  al., 
1987b].  The  SnapCon  graphics  package  is  described  in  Section  3.3. 

3.1  Input  Data 

(1)  Control  Data 

The  control  data  determine  what  type  of  calculations  are  to  be  performed.  There  is 
only  one  input  data  file  in  this  category  (the  control  file,  Figure  3.1).  NetSim  prompts 
for  the  name  of  the  control  file  at  the  start  of  each  run.  The  user  is  given  an  opportun¬ 
ity  to  view  and/or  modify  the  contents  of  this  file,  write  any  changes  to  disk,  and  to 
proceed  by  reading  the  rest  of  the  input  data  files. 

The  control  data  include  the  NetSim  run  type  which  is  either  detection ,  location ,  or 
detloc  (detection  and  location).  If  the  run  type  is  set  to  detection,  then  the  location 
module  is  not  activated.  However,  if  the  run  type  is  set  to  location ,  the  detection 
module  must  calculate  the  individual  station  probabilities  since  these  are  needed  to 
estimate  location  uncertainty.  If  the  run  type  is  set  to  detloc  then  both  modules  are 
activated,  and  the  location  uncertainties  for  each  epicenter  are  calculated  at  the 
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Figure  3.1.  Schematic  top-level  data  flow  diagram  for  NetSim. 


detection  threshold  for  the  network.  The  run  subtype  is  also  included  in  the  control 
input.  This  is  either  threshold  (to  calculate  the  detection  threshold  at  fixed  confidence 
level),  or  probability  (to  calculate  the  probability  of  detection  for  fixed  event  size). 
Note  that  event  size  can  be  represented  in  terms  of  explosion  yield,  seismic  moment, 
mb,  Ms,  or  MLg.  The  conversions  between  these  measures  are  included  in  the  input 
data  for  each  source  medium.  Other  control  data  include  the  detection  criteria  (e.g., 
the  wave  types  and  number  of  stations  required  for  detection),  probability  of  detection 
for  the  network  (for  a  threshold  run),  event  size  (for  a  probability  or  location  run), 
location  confidence  level,  number  of  Monte  Carlo  realizations  (=1  for  probability¬ 
weighting),  frequencies  to  inc’*ide  in  the  calculation,  and  the  names  of  the  input  data 
files. 

(2)  Source  Data 

The  source  data  are  used  to  calculate  the  excitation  of  each  wave  as  a  function  of  fre¬ 
quency  and  source  medium,  Sjtf).  There  are  four  input  file  types  in  this  category:  epi¬ 
center  grid  file,  source  medium  file,  source  medium  grid  file,  and  the  source  spectra 
files. 

The  epicenter  grid  file  specifies  the  event  epicenters  at  which  to  calculate  the  detection 
thresholds,  probability  of  detection,  or  location  uncertainties.  The  source  medium  file 
specifies  source  scaling  and  excitation  factors  for  each  source  medium.  This  includes 
the  medium  name  (e.g.,  granite,  tuff,  salt),  the  file  name  for  the  source  spectrum,  the 
density  and  velocity  of  the  source  medium,  the  scaling  relations  between  log  moment 
and  other  measures  of  event  size  (yield,  mb,  M„  and  MLg),  and  the  excitation  factors 
for  each  wave,  k*.  The  source  grid  file  identifies  the  geographic  locations  of  each 
source  medium.  The  source  spectra  for  each  source  medium  are  stored  in  two- 
dimensional  data  files  that  tabulate  M0/f)  as  a  function  of  scalar  seismic  moment  and 
frequency. 

(3)  Propagation  Data 

The  propagation  data  are  used  to  determine  the  frequency-dependent  attenuation  B^Aij) 
and  travel  times  (and  their  standard  deviations)  for  each  wave  as  a  function  of  path 
medium.  There  are  four  input  file  types  in  this  category:  path  medium  file,  path 
medium  grid  file,  attenuation  files,  and  travel  time  files. 

The  path  medium  file  specifies  attenuation  and  travel  time  parameters  for  each  path 
medium.  This  includes  the  medium  name  (e.g.,  tectonic  or  stable),  file  names  for 
attenuation  and  travel  times  for  each  wave,  the  standard  deviation  of  the  attenuation 
curves  for  each  wave,  the  travel  time  and  azimuth  standard  deviations,  and  the  source 
and  receiver  amplitude  correction  factors  for  teleseismic  distances.  The  path  grid  file 
identifies  the  geographic  locations  of  each  path  medium.  The  attenuation  3^)  for 
each  wave  and  path  medium  is  tabulated  in  two-dimensional  data  files  as  a  function  of 
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distance  and  frequency.  The  travel  times  are  tabulated  as  a  function  of  depth  and  dis¬ 
tance. 

(4)  Site/Station  Data 

The  site/station  data  include  the  local  site  response  for  each  station  SR&  (which 
includes  the  free  surface  amplification,  instrument  response,  and  the  frequency- 
dependent  beam  gain  for  array  stations),  and  individual  station  parameters.  There  are 
two  input  files  in  this  category:  station  file  and  site  response  file. 

The  station  file  includes  the  names  and  geographical  coordinates  of  each  station,  file 
names  for  ambient  noise  spectra,  the  standard  deviation  of  log  noise,  station  reliability, 
amplitude  station  corrections,  signal-to-noise  ratio  required  for  detection  of  each  wave, 
station-specific  standard  deviations  (travel  time,  azimuth,  and  log  amplitude),  file 
names  for  coda  decay  rates,  and  file  names  for  the  site  response  for  each  wave  and  sta¬ 
tion.  The  site  response  is  tabulated  in  a  one-dimer.sional  data  file  as  a  function  of  fre¬ 
quency. 

(5)  Noise  Data 

The  noise  data  define  the  ambient  station  noise  spectra  and  the  coda  decay  rates  used 
to  calculate  noise  spectra  for  secondary  phases.  There  are  two  input  files  in  this 
category:  the  ambient  noise  file  <c  the  noise  factor  file  for  secondary  phases.  Both  of 
these  files  are  indexed  in  the  jcatioi,  file. 

There  is  a  separate  ambient  noise  file  for  each  station.  These  files  tabulate  ambient 
noise  power  spectral  density  as  a  function  of  frequency.  The  noise  factor  file  for 
secondary  phases  tabulates  the  coda  decay  rate  of  the  previous  arrival  as  a  function  of 
distance,  y*(Ay). 

3.2  Output  Data 

NetSim  writes  two  types  of  output  files;  a  graphics  file  used  by  SnapCon  to  display  the 
results  as  contours  on  a  map  (Section  3.3),  and  a  list  file  that  contains  detailed  infor¬ 
mation  about  the  detection  and  location  parameters  for  each  epicenter  and  for  each  sta¬ 
tion.  The  graphics  file  is  called  prefix.gra  and  the  list  file  is  called  prefix.lst ,  where 
prefix  is  specified  in  the  control  input 

The  graphics  file  includes  a  table  that  lists  event  size,  probability  of  detection  for  the 
network,  and  location  uncertainties  for  each  epicenter  in  the  grid.  For  threshold  runs, 
the  event  size  is  equal  to  the  detection  threshold  for  the  network.  Otherwise,  it  is 
equal  to  the  fixed  value  specified  in  the  control  input  (e.g.,  for  probability  or  location 
runs).  The  location  uncertainties  are  reported  as  the  length  of  the  semi-major  and 
semi-minor  axes  of  the  epicenter  location  error  ellipse,  the  strike  of  the  semi-major 
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axis,  and  the  depth  uncertainty.  For  Monte  Carlo  simulations,  the  fraction  of  simula¬ 
tions  with  enough  data  to  locate  the  event  is  also  reported  in  the  graphics  file. 


The  list  file  includes  all  of  the  information  contained  in  the  graphics  file  and  informa¬ 
tion  about  the  detection  and  location  capabilities  of  individual  stations  in  the  network. 
The  list  file  includes  the  following  information  for  each  station,  wave  and  epicenter: 
detection  probability  P^  epicentral  distance  and  event  azimuth,  frequency  of  the  max¬ 
imum  signa\-to-noise  ratio,  standard  deviation  of  log  amplitude,  azimuth  and  arrival 
time  standard  deviation,  and  depth  and  epicenter  importances  [Ciervo  et  al. ,  1985], 

3.3  Graphics 

SnapCon  is  a  graphics  package  developed  by  3 ran  et  al.  [1987b]  to  plot  SNAPJD 
results  as  contours  on  regional  and  world  maps.  Normally,  SnapCon  is  used  to  con¬ 
tour  detection  thresholds  or  location  uncertainties.  However,  it  can  also  be  used  to 
contour  the  difference  or  percentage  change  between  two  simulation  runs.  Using  these 
options,  it  is  possible  to  estimate  the  change  in  detection  capability  due  to  changes  in 
input  parameters  such  as  network  geometry,  detection  criteria,  or  attenuation  curves. 
SnapCon  uses  the  NCAR  (National  Center  for  Atmospheric  Research)  graphics  library 
with  enhancements  by  MINEsoft,  Ltd.  SnapCon  can  be  used  to  generate  one  or  more 
plots  during  a  single  run,  and  these  plots  are  stored  in  a  metacode  graphics  output  file. 
These  metacodes  can  be  plotted  on  selected  output  devices  (e.g.,  SunView  or 
PostScript)  using  metacode  translation  programs. 

The  NetSim  output  graphics  file  has  a  slightly  different  format  than  the  SNAP/D  output 
graphics  file,  so  we  made  minor  modifications  to  SnapCon  so  that  we  could  use  it  to 
contour  NetSim  results.  This  new  version  of  SnapCon  has  the  same  functionality  as 
the  version  described  in  the  SnapCon  user’s  and  programmer’s  guide  written  by  Bratt 
et  al.  [1987b].  Examples  of  SnapCon  graphics  are  given  in  Section  4  of  this  report. 

Bratt  et  al.  [1987c]  describe  a  graphics  package  (also  based  on  the  NCAR  graphics 
library)  to  plot  SNAP/D  results  stored  in  the  list  file.  This  package,  called  PRIMP, 
plots  the  probability  of  detecting  each  wave  at  individual  stations  for  a  given  epicenter 
(Pyt),  and  location  importances  which  measure  the  relative  contribution  of  a  particular 
datum  to  event  location.  These  plots  are  useful  for  determining  detection  and  location 
capabilities  of  individual  stations  in  the  network.  However,  we  have  not  yet  modified 
the  PRIMP  graphics  package  to  read  the  list  file  generated  by  NetSim. 
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4.  SAMPLE  RUNS 


This  section  presents  sample  NetSim  and  SnapCon  runs  for  frequency-independent  and 
frequency-dependent  input.  The  simulations  presented  in  this  section  are  only  illustra¬ 
tive  examples  of  the  NetSim  computer  program.  These  simulations  are  not  intended  to 
represent  accurate  estimates  of  the  detection  and  location  capability  of  seismic  net¬ 
works  in  the  Soviet  Union.  Instead,  the  examples  in  this  section  illustrate  the  main 
functionality  of  NetSim,  the  importance  of  frequency  dependence  for  simulating  capa¬ 
bility  at  regional  distances,  and  the  inability  of  the  probability-weighted  ( SNAP/D ) 
approach  to  properly  account  for  the  effect  of  undetected  phases  on  location  uncertain¬ 
ties  for  events  near  the  detection  threshold  of  the  network. 

Section  4.1  presents  sample  NetSim  runs  at  fixed  frequency.  The  results  for  two 
separate  calculations  (1  Hz  and  10  Hz)  are  compared.  Included  in  that  section  are  list¬ 
ings  of  the  input  and  output  data  files,  and  a  complete  transcription  of  an  interactive 
NetSim  session.  Section  4.2  presents  simulations  for  frequency-dependent  input,  and 
compares  them  to  the  simulations  for  fixed  frequency  presented  in  Section  4.1.  Sec¬ 
tion  4.3  compares  location  uncertainties  computed  using  the  probability-weighted 
approach  to  those  computed  using  the  Monte  Carlo  approach. 

4.1  Frequency-Independent  Input 

This  section  describes  sample  NetSim  calculations  for  fixed  frequency.  These  calcula¬ 
tions  have  been  compared  to  SNAP/D  calculations  to  verify  that  these  programs  give 
the  same  output  when  they  are  given  identical  (frequency-independent)  input.  The 
main  purpose  of  this  section  is  to  give  examples  of  the  input  and  output  files  for  Net¬ 
Sim,  and  to  illustrate  the  NetSim  data  handling  and  user-interface  modules. 

INPUT  DATA 

The  input  parameters  described  here  are  frequency-dependent.  However,  the  simula¬ 
tions  presented  in  this  section  are  calculated  for  a  single  frequency,  so  they  are  the 
same  as  simulations  for  frequency-independent  input  (e.g.,  like  SNAPID  simulations). 
The  advantage  of  using  frequency-dependent  input  data  files  is  that  they  can  be  used 
without  modification  to  simulate  broad  band  detection  capability  (as  demonstrated  in 
the  next  section)  or  detection  capability  at  any  fixed  frequency.  The  sample  input  data 
files  in  each  of  the  five  categories  (control,  source,  propagation,  site/station,  and  noise) 
are  described  below.  None  of  the  sample  runs  include  secondary  phases,  so  the  only 
wave-dependent  parameters  included  in  the  sample  input  files  are  for  P  phases.  These 
sample  data  files  are  based  on  our  experience  with  the  NORESS  array  in  Norway 
[Sereno  et  al.,  1988].  However,  we  have  extrapolated  our  NORESS  results  to  much 
larger  distances  and  somewhat  higher  frequencies.  Therefore,  the  simulations  have 
large  uncertainty  (particularly  for  events  near  the  periphery  of  the  network)  and  should 
not  be  interpreted  as  accurate  estimates  of  detection  capability  in  the  Soviet  Union. 
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Control  Data.  Table  4.1  lists  the  contents  of  the  sample  control  file.  This  file 
specifies  that  the  90%  MLg  thresholds  for  detecting  3  P  phases  at  1.0  Hz  will  be  calcu¬ 
lated.  Magnitude  search  bounds  are  1.0  and  8.0.  The  names  of  other  input  files  are 
listed  at  the  bottom  of  the  control  file. 

Source  Data.  Table  4.2  lists  the  contents  of  the  four  source  data  files  (epicenter  grid, 
source  media,  source  grid,  and  source  spectra).  These  files  specify: 

•  A  15°  x  15  0  epicenter  grid  for  Eurasia 

•  Homogeneous  granite  source  medium  grid 

•  Source  medium  properties  (granite) 

Density  =  2500  kg/m3 

-  Compressional  wave  velocity  =  5000  m/s 

•  Source  scaling  relations  (granite) 

-  log  Mq  =  1.08  MLg  +  10.6 

-  P-wave  excitation  factor,  k  =  1.0 

-  Source  spectra  as  a  function  of  log  M0  (Figure  4.1) 

Propagation  Data.  Table  4.3  lists  the  contents  of  the  four  propagation  data  files  (path 
media,  path  grid,  P-wave  attenuation,  and  P- wave  travel  times).  These  files  specify: 

•  Homogeneous  propagation  characteristics  (stable) 

•  Window  length  for  P  wave  spectra  =  5.0  s 

•  Standard  deviations  for  P  phases 

-  Log  amplitude  standard  deviation  =  0.3 

-  Arrival  time  standard  deviation  =  1.5  s 

-  Azimuth  standard  deviation  =  8.0° 

•  P-wave  attenuation  versus  frequency  (Figure  4.2) 

•  P-wave  travel  times 

Site/Station  Data.  Table  4.4  lists  the  contents  of  the  two  site/station  data  files  (sta¬ 
tion,  site  response).  These  specify: 

•  A  network  of  20  internal  and  13  external  arrays  (Figure  4.3) 

•  The  same  ambient  noise  spectra  and  standard  deviation  for  all  stations 

•  A  P-wave  SNR  for  detection  =  4.0  for  all  stations 

•  Receiver  medium  properties  for  all  stations 

Density  =  2500  kg/m3 
Compressional  wave  velocity  =  5000  m/s 
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CONTROL  DATA 


CONTROL  FILE:  ./flrunl.ent 


Frequency-Independent  Detection  Thresholds 
Flrunl 
. /results 

•  / 

detection 

threshold 

(P/3) 

P 

mlg 

3.000  1.000  8.000  0.900 

0.900  100  0.500  0.900 

fixed  0.000 
1 

1.000 
euras.epi 
granite,  sad 
granite. sgr 
stable,  pad 
stable. pgr 
a20xl3.sta 


0.000 


20.000 


Table  4.1.  Sample  control  data 
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SOURCE  DATA 


EPICENTER  PUB:  ./euraa.epi 

Eurasia  Epicenters 
7  13 

0.000  0.000 

15.000  15.000 


SOURCE  MEDIA  FILE:  . /granite . sad 

Granite  Source  Madia 
. /arctab 
1 

granita 

expl.aor 


2500.000 

0.000 

1.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

1.080 

0.000 

10.600 

0.000 

5000.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

o.coo 

SOURCE  GRID  FILE:  . /granita. agr 


Granita  Sourca  Madia  Grid  (5x5  dagraa) 


36 

-90.000000 

5.000000 

granita 

granita 

granita 


72 

-180.000000 

5.000000 

-90.000000 

-90.000000 

-90.000000 


-180.000000 

-175.000000 

-170.000000 


granita 

granite 

granite 


S5.r.o000O 

85.000000 

85.000000 


165.000000 

170.000000 

175.000000 


SOURCE  spectra  FILE:  . /arctab/axpl . aor 


Explosion  Source 
20 

1.000  2.000  3.000 

9.000  10.000  11.000 

17.000  18.000  19.000 

6 

11.000  12.000  13.000 

Sourca  amplitudes  for  frequency 
11.000  12.000  13.000 

Sourca  amplitudes  for  frequency 
11.000  12.000  13.001 


4.000  5.000 

12.000  13.000 

20.000 

14.000  15.000 

1.00000 

14.001  15.003 

2.00000 

14.002  15.002 


Sourca  amplitudes  for  frequency  «  18.0000 

11.002  12.003  12.910  13.401  13.737 

Sourca  amplitudes  for  frequency  *  19.0000 

11.002  12.002  12.890  13.356  13.690 

Sourca  amplitudes  for  frequency  m  20.0000 

11.002  12.002  12.867  13.312  13.645 


6.000 

14.000 

20.000 

17.912 

17.310 

15.401 

15.354 

15.310 


7.000  8.000 

15.000  16.000 


Table  4.2.  Sample  source  data. 
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Figure  4.1.  Input  source  spectra  as  a  function  of  scalar  seismic  moment,  Af0.  Each 
curve  is  labeled  by  log  Mq.  The  comer  frequency  scales  inversely  with  the  cube  root 
of  the  long-period  level.  The  source  scaling  is  from  Sereno  et  al.  [1988].  These 
spectra  are  stored  in  a  two-dimensional  data  table  called  expUor  (Table  4.2). 


PROPAGATION  DATA 


PASS  MEDIA  FILE:  ./atabla. pod 


Stabla  Path  Madia 
. /path tab 


5.000 

-1.000 

-1.000 

0.000 

-1.000 

-1.000 

0.000 

-1.000 

-1.000 

0.000 

-1.000 

-1.000 

0.000 

-1.000 

-1.000 

0.000 

-1.000 

-1.000 

0.000 

-1.000 

-1.000 

0.000 

-1.000 

-1.000 

0.000 

-1.000 

-1.000 

o.coo 

-1.000 

-1.000 

2 

rafaranca 

patabla. atn 

patabla . tin 

0.300 

0.000 

0.000 

1.500 

8.000 

“  0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

~  0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

~  0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo” 

0.000 

0.000 

atabla 

patabla.atn 

patabla. tin 

0.300 

0.000 

0.000 

1.500 

8.000 

0.000 

0.000 

0 . ooo- 

0.000 

0.000 

“  0.000 

0.000 

0 . ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

0 . 000~ 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

0.000 

0.000 

o.ooo- 

0.000 

0.000 

“  0.000 

0.000 

o.ooo- 

0.000 

0.000 

PATH  GRID  PILE:  . /atabla. p?r 

Stabla  Path  Madiua  Grid 
36  72 

-90.000000  -180.000000 

5.000000  5.000000 


atabla 

-90.000000 

-180.000000 

atabla 

-90.000000 

-175.000000 

atabla 

-90.000000 

-170.000000 

atabla 

85.000000 

165.000000 

atabla 

85.000000 

170.000000 

atabla 

85.000000 

175.000000 

Table  4.3.  Sample  propagation  data. 
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insmOXZON  FILE:  . /patht.vb/p»table . atn 


stable  -  p 
20 


0.500 

1.500 

2.500 

3.. 500 

4.500 

5.500 

6.500 

7.500 

8.500 

18.500 

50 

0.000 

9.500 

17.500 

10.500 

18.500 

11.500 

20.000 

12.500 

13.500 

14.500 

15.500 

2.000 

4.000 

6.000 

8.000 

10.000 

12.000 

14.000 

16.000 

18.000 

20.000 

22.000 

24.000 

26.000 

28.000 

30.000 

32.000 

34.000 

36.000 

38.000 

40.000 

42.000 

44.000 

46.000 

48.000 

50.000 

52.000 

54.000 

56.000 

58.000 

60.000 

62.000 

64.000 

66.000 

68.000 

70.000 

72.000 

74.000 

76.000 

78.000 

80.000 

96.000 

Attenuation 

82.000  84.000 

98.000 

for  frequency  0 . 

86.000 

.500000 

£8.000 

90.000 

92.000 

94.000 

0.000 

-6.138 

-6.617 

-6.933 

-7.182 

-7.395 

-7.585 

-7.759 

-7.922 

-8.076 

-8.222 

-8.363 

-8.499 

-8.631 

-8.760 

-8.886 

-9.010 

-9.131 

-9.251 

-9.368 

-9.484 

-9.599 

-9.712 

-9.824 

-9.936 

-10.046 

-10.155 

-10.263 

-10.371 

-10.478 

-10.584 

-10.690 

-10.795 

-10.899 

-11.003 

-11.106 

-11.209 

-11.312 

-11.414 

-11.516 

-11.617 

-11.718 

-11.819 

-11.919 

-12.019 

-12.119 

-12.219 

-12.318 

-12.417 

-12.516 

. 

• 

. 

. 

. 

. 

Attenuation 

for  frequency  20.0000 

: 

: 

: 

: 

0.000 

-6.644 

-7.628 

-8.450 

-9.206 

-9.925 

-10.620 

-11.300 

-11.969 

-12.628 

-13.281 

-13.927 

-14.569 

-15.208 

•15.842 

-16.474 

-17.104 

-17.731 

-18.356 

-18.979 

-19.601 

-20.222 

-20.841 

-21.459 

-22.076 

-22.692 

-23.307 

-23.921 

-24.535 

-25.147 

-25.760 

-26.371 

-26.982 

-27.592 

-28.202 

-28.811 

-29.420 

-30.028 

-30.636 

-31.244 

-31.851 

-36.698 

-32.458 

-37.302 

-33.065 

-33.671 

-34.277 

-34.882 

-35.488 

-36.093 

TRAVEL 

TIME  TABLE:  . /patbtab/pstable . tin 

Travel-tine 

9 

v*  distance  and  depth  -  Pn  stable 

0.000 

15.000 

181 

0.000 

1.000 

2.000 

3.0C0 

4.000 

5.000 

6.000 

7.000 

8.000 

9.000 

10.000 

11.000 

12.000 

13.000 

14.000 

15.000 

168.000 

169.000 

170.000 

171.000 

172.000 

173.000 

174.000 

175.000 

176.000 

177.000 

178.000 

179.000 

180.000 

Travel-tine 

for  source  depth 

0. 

0.000 

21.134 

34.863 

48.582 

62.291 

75.988 

89.670 

103.335 

116.978 

130.597 

144.189 

157.752 

171.289 

184.775 

198.199 

211.480 

1206.987 

1207.651 

1208.261 

1208.815 

1209.314 

1209.755 

1210.140 

1210.467 

1210.736 

1210.946 

1211.158 

1211.253 

1211.288 

Travel -tine 

for  source  depth 

15.0000 

2.500 

17.585 

33.185 

46.902 

60.610 

74.305 

87.985 

101.646 

115.286 

128.902 

142.490 

156.048 

169.579 

183.058 

196.466 

209.722 

1204.488 

1205.153 

1205.762 

1206.317 

1206.815 

1207.257 

1207.642 

1207.969 

1208.238 

1208.448 

1208.621 

1208.713 

1208.745 

Table  4.3.  Sample  propagation  data  (continued). 
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P-wave  Attenuation  Curves 


4  8  12  16  20 

Epicenteral  Distance  (Degrees) 


Figure  4.2.  Input  P- w»5ve  attenuation  curves  as  a  function  of  frequency.  These  curves 
are  based  on  attenuation  along  paths  to  the  NORESS  array  [Sereno  et  al.,  1988]. 
These  B(A,f)  are  stored  in  a  two-dimensional  data  table  called  pstable.atn  (Table  4.3). 
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SITE /STATION  DATA 


STATION  FILE:  ./a2 0x13. at a 


13  X  AR,  20  IN  AR 
.  /stf.tion 
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HOR  ambient . noi 


60.735 

11.542 

2500.000 

0.000 

0.000 

O.COO 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

SHIO  ambient. noi 

25.567 

91.883 

2500.000 

0.000 

0.000 

0.000 

0.000 

0.000 

o.occ 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.100 

1.000 


4.000 

0.000 

0.000 

5000.000 

array. ara 

0.000 

0.000 

0.000 

C.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

““ 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

O.COO 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.^00 

0.100 

1.000 

4.000 

0.000 

0.000 

5000.000 

array. arc 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

o.ooo 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.  Q00 

0.000 

0.000 

0.000 

21  ambient. noi 

0.100 

77.500 

105.100 

1.000 

2500.000 

0.000 

0.000 

4.000 

0.000 

0.000 

5000.000 

array. ar a 

0.000 

0.000 

0.000 

o.ooo 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

*■* 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

7.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

- 

SITE 

RESPONSE  FILS 

>:  ./station/array 

.srs 

to  Reapcnse 

,  array  station 

11 

1.000 

0.699 

3.000 

0.661 

5.000 

0.619 

7. 090 

0.573 

9.000 

0.521 

11.000 

0.462 

13.000 

0.394 

15.000 

0.314 

17.000 

0.215 

19.000 

0.086 

20.000 

0.004 

Table  4.4.  Sample  site/station  data. 
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60.0° 


SIMULATION  NETWORK 

-40.0®  20.0° 


80.0° 


oFrd«e^r^  ~  used 

the  station  file,  a20xl3.sta  (Table  4.4).  tUl  coordinates  we  stored 


•  Site  response  spectrum  for  all  stations  (Figure  4.4) 

Noise  Data.  Table  4.5  lists  the  contents  of  one  of  the  two  noise  data  files  (ambient 
noise).  The  other  noise  input  file  (noise  factor  for  secondary  phases)  is  not  needed  for 
this  sample  calculation  since  the  detection  criteria  only  include  P  waves.  The  ambient 
noise  file  is  the  power  spectral  density  for  all  stations  in  the  network  (Figure  4.5). 

NETSIM  TRANSCRIPTION 

A  complete  transcription  of  an  interactive  NetSim  session  that  uses  the  input  data  files 
described  above  starts  on  page  31  of  this  report.  Comments  are  separated  from  the 
transcription  by  a  vertical  line  on  each  page.  This  transcription  is  included  as  exam¬ 
ples  of  the  NetSim  data  handling  and  user-interface  modules. 
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Local  Site  Response  (Array  Stations) 


Figure  4.4.  Input  local  site  response  for  array  stations.  Only  the  aiTay  gain  portion  of 
the  site  response  is  included.  We  assume  a  linear  decrease  in  beam  gain  from  VJv  at  1 
Hz  to  1.0  at  20  Hz  [Screno  et  al.,  1988].  The  site  response  ?.s  a  function  of  frequency 
is  stored  in  a  one-dimensional  data  table  called  array. srs  (Table  4.4). 
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NOISE  DATA 


AMBIENT  NOISE  FILE:  . /»t»tion/«mbi«nt .aoi 

Ambient  Noi*a,  Noras* 

253 

0.234  5.530 

0.312  4.410 

0.391  3.620 


19.766  -5.550 
19.844  -5.520 
19.922  -5.460 


Table  4.5.  Sample  noise  data. 
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Ambient  Noise  -  NORESS 


Figure  4.5.  Input  ambient  noise  power  spectral  density.  The  same  ambient  noise 
spectrum  is  assumed  for  all  stations  in  the  network.  This  is  the  average  single-channel 
noise  spectrum  at  NORESS  [Suteau-Henson  and  Bache,  1988].  The  ambient  noise  is 
stored  in  a  one-dimensional  data  table  called  ambient. noi  (Table  4.5). 


30 


II 

AJ  C 

m  m 

*H  ®  . 

H  XJ  = 
AJ  » 
A>  C  ® 

s.2  r- 

2  3.  tn 


£  5 

°  -  C 

®  H 

X2  9 

C  -H 

X3  -H 

AJ  « 

•  -H  M 

AJ 

*rl 

X3  • 

14  «H 

A>  *0  AJ 

0  = 

•  C 

•w  o 

«H 

tn  aj  • 

*H  i4 

« 

C  ® 

^  AJ 

-rf  o  N 

S3-0 

«  C  JJ 

H  *rt  H 


«  • 

or 

no' 

of- 

M 

H 

•  p 

-H 

<H  :  *0 

ss 

>1  j§ 

4  (4 
♦H 

do® 
•  AJ  X- 
>  AJ 

« 

o 

ft 

Ojaj  m 

u 

«  • 

tJ  ® 

*0  H 

f-i  • 

•  P  V 
■H  4  i4 

o  c  • 
e  o  « 

.  0  • 

M  r~f 

O  Jt  ^ 

•  *4 

*  c 

«  4 

«H  4  A> 

>  4 

®  •  *0 

H  AJ 
■H  M 

*  M  « 
0  *  XJ 
C  AJ 

«.  o 
>1  f 
4  H  ® 

g  X3 

«W  AJ 
14 

•  H  Ol 
HOC 
P  H 

a>  *0 

•  G  * 

x:  o  « 

P  O  M 


.  4  * 

5  HO. 

4  «*4  4 

5  *43 

•  4  . 

rH  4J  U 

*H  4  0 

*M  *0  ^ 

H  M  4 

0  C  AJ 

14  X3  4 

AJ  AJ  *0 

C  O 

O  • 

O  14  r-4 

O  *H 
>t  »M 

•H  «  H 

•o  ®  o 

O  B  14 

g  2  AJ 

^  c  c 

•  o 

•  •  V 


C  4J 

®  4 

&  S 

®  ft  0 

u  ®  4 

*4  4  0 

®  ®  & 

gss. 


o  C  w 
M  9  4 
a>  o  ft 

o  ®  ® 

o  «  x: 

®  aj 
•  X3 
XS  AJ  *M 

AJ  0 


AJ 
<v.  C 
AJ  W  U 


fa  c 

C  H  B 

s  ss 


•  •  <n 

4J  M\ 


o  o  o  o 
o  o  o  o 
o  o  o  o 
o  o  o  o 
o  o  ©  o 
•  •  •  o» 

n  H  CD  . 

o 


J  XJ  ft.  H 

■d  aj  ~  *  a 


o  o  o  o 

o  o  o  o 

o  o  o  o 

©  o  o  o 

•  o  o  o 

o  in  o\ 
fH  •  O  •  •  TJ 
o  o  o  o  •  • 

O  H  X  O  « 


H  Q 

S  2  a 

K 


•  •  o 

:  s  s  * 

c  fifli 

I  W  -r< 

AJ  AJ  K 

!  C  C  M 
I  O  O  K 

:  U  U  a 


4ffiUPMX 


q  ®  4  as 

«  3^*5 

a  fc&i* 

O  •  55  § 
r*  HP 
O  AJ  U*  •  H 
U  fl  ®  H  H 
AJ  O  H  «H  K 
C  U  In  Pm  ^ 
o 

U  .... 
*  u  X 


•  u  u  i  jis  >  I;  •  • 

S.33ff  85S3SS2 

sacpfl«o*<>s>M 

oouuat>H«*««B< 


„  0  XI 

AJ  AJ  0  H  <H  I 

•  I  O  M  0HW1 
V  0  •  r-4  4  H  X3  X3  4 
X4  H  *0X1  0  AJ  AJ  AJ  1 

04J  |  o  o  o«i 

M  3  *  U  P  H  M  • 
fa  0  M  0.  C  0.  fa*0  < 


KnuoMhonHbx^xxooiOi^cQ 


a 

x  US 

a  °5 

a  m  a 
•  ^  ^  « 

^  8S* 

•A  M  5 
0  AJ  &*• 
H  6  •  H 
AJ  O  M  *H 
C  U  N  Cm 
o 

u  ... 

*<  «  u 


EXIT  TO  MAIN  MEND 


ft 

• 

P 

H 

? 

S' 

c 

I 

u 

%4 

3 

i 

* 

ft 


43 

41 


c 

o 

■H 

41 

ft 


ft 

o  • 

•  5 

JO  03 
41 


■d 

c  « 

ft  4J 
3  « 

• 

2  & 


C 

o 

•H 

41 

ft 

0 


4J 

P 

O 

43 

41 


3 

C 

•H 

4> 

C 

O 


43 

P 

& 

-H 

O 


c 

H 

>1 


9  U 

O 

O 

> 

9 

> 

« 

•H  9 

0 

> 

ft 

M 

V 

41 

9 

C 

U 

0 

H 

A 

04 

0 

4> 

S’ 

ft 

41 

t 

9  9 

o 

43 

: 

A 

a  g. 

V 

41 

V 

• 

O 

9 

9  U 

« 

41 

ft 

43 

A  *M 

M 

41 

0 

41 

9 

• 

0 

9 

• 

41 

E 

ft 

r-K 

«  c 

c 

r4 

9 

41  O 

9 

3 

ft 

9 

« 

41 

ft 

•H  >i 

9 

ft 

u 

H  H 

U 

u 

U 

9 

^  c 

9 

8 

ft 

9 

6  o 

9 

ft 

P 

•H 

3 

•H 

3 

CO  c 

SO 

9 

41  • 

« 

41 

ft 

A 

9  9 

A 

ft 

F* 

2  0 

H 

2 

43 

l 

c 


4> 

c 

ft 

M 

H 


ft 

43 

43 

m 

43 


to 

43 

9 

as 


c 

o 

V 

ft 

a 


•p 

o 

B 

43 

P 

o 

JO 

43 


g 

1  5 

ft  p 

P  c 

l 

6  VM 

5  g 

C 

o 

o 

o 

43 


43 

c 

ft 

m 

u 

9 


U 

su 

9 

JO 

43 

O 

43 

9 

C 

u 

3 

43 

2 


CO 

43 


41 

p 

H 

c 

3 

i 

0 

ft 

V 

■H 

ft 

«H 

43 

41 

H 

O 

41 

(4 

0< 

41 

ft 

C 

u 

o 

u 

0 

9 

c 

0 

•H 

41 

A 

ft 

P3 

43 

D 

41 

V 

0 

41 

m 

u 

9 

43 

c 

9 


U 

S  8 

3  n 

•  43 


E 

P 

43 

9 

U 


43 

I 


u 

9 

A 

43 

O 

9 

A 

43 

&» 

c 

* 

n 

g 

■H 

43 

C 

O 

a 

o 

43 


43 

V 

9 

H 

S. 

uS 


3- 


A 

43 


43  «H 

C 

• 

M  ft 

M  V 

B-S 

•  4J 

ja  a 
u  c< 
c 

•  M 

43 

c  c 

T!fl 

N 

6  o 

«rt 

*3 

2  C 


C 

*  o 

ft  *ri 

9  43 
f-4  1 4 

•H  U 

%4  *H 
*  +4 


S 

ft  ft 

r-l  41  ■© 
•H  ft  O 

& 

•M  *0  8 

O 

H  41  H 

41 

0  3ft 

O 

0  0441 

ft 

4>  C  *M 

M 

C  <H  ft 

5 

O 

0  ft  * 

43  • 

ft 

%  41  O' 

H 

ft 

•ri 

0  <H  ft 

o 

ft 

H 

>1  « 

O 

43  C  ft 

^  «M<U 

o 

U  9 
9  m  43 

>  43  9 

O  C  *0 
a 

4l  41  41 
0  C  P 

k  o  a« 

D  U  £ 
9  -H 

ft 

*  43  • 
q  4i  A 

C  ** 

as® 

u  e  c 

©  ^  9 

ft  » 

3  •  • 

«H  43 
>  *H 

9  M 

H  » 


I 


O 

O 

o 

o 

o 


•I 

• 

3 
H 

4 

> 

S’ 

C 


2 


2 

H 

H 

§ 

O 

o 

41 


M 

O 

I 

g 

c 

43 

9 

43 


g 

H 

<0 


ft 

9 

Q 

• 

H 


O 

M 

41 

a 

o 

u 


9 

43  9 

ft  • 

a  P 

*H  « 

sn. 

O  9  2  K 
J  &  •  n 

§  S3 s 

UhkM 


*  a  u  k 


p 

SB  O 

!  B 


2 

H 

H 

§ 

U 

o 

41 

I 


V  u  u 

8  ©'•d  u  9  o 

«H  ft  •  S  0*41 

S’.-  iWm.  8. 

.  43  43  9  9  <3  C 

•  ft 

1  2  S  S  -3  -a  s  -s 

a  3  R  *4  43  43  04 

C  ft  0*  0*  ft  ft  ft  o 

43 

5  .  e 

*  frlmg  5 

H  '••••• 

M  HH  WHHH  t) 

4)  «H  «H  <H  *H  «H  ft 

C  %4  <M  <M  «M  %4  <M  r-4 

0  • 

u  ••••••  CO 

<  m  V  Q  M  to 


b 


9 

a 

p 


9  • 


0 

N 

C 

O 

u 


u 


Afflux 


- 

p 


a 


i 

§ 

u 


H 

P 

i 

& 


•  9 
9  « 

H  • 

•H  *H 

#  -I 

•H  4  « 

jJ  °  1g 

•  IMP 

2  8  2  ^  & 

H  CH 

O  O  O  t*4> 

n  u  n  C  P 
4)  4>3  o 

§&gf5: 

‘  lu«Si 

•  x  •  _ _ 

2  2  g 

,  >  C  C  ft 
M  #  «  «H  *4 
9  M  43  43  43  K 
4i  >  -H  C  C  M 
C  •  M  O  O  X 

mmuum 


X  «  u  Q  W  X 


s 


?  d»*d  u  9 

•H  •  «  5  &  43 

«  «  •  T  T  • 

.  43  43  9  9  <n 

.  2  S  S-S-35 

3  M  H  43  43  C4 

9  D»  t*  •  •  • 


i 


(4 

u 

a 

s 


•HH  • 

•  ft  «H  H  • 

n  ft  3  H 

•H  <H  • 

ft  3  *0  ^  ft  H 

sO  *H  4  «rt 
^  H  -H  TJ  Pm 

•  2  e  *0  «H  . 

4i  •  n  a 
G  9  9  X  O  O 

f  CJ  O  -H 

OMM4J4J 
«H  P  3  41  44  9 
0U  0  O  4  4  41 
ft|  CO  CO  {ft  fW  <0 

^  ^ 
fi  C  C  C  fi  fi 

5  t1  *3  *3  3  ?} 

9  9  9  9  9  9 

2  2  5  2  2  2 


o 

M 

41 

O 

u 

ft 

• 

2 

M 

• 

41 


32 


TRANSLATING  DETECTION  CRT 


-  u 

4  9 

W  C 

*  8 

4  «rl  -H 

«  0  S* 

^  C  JG 
41  O  O 
P  «H  4 

•H  CM 
44  0 
J3  «  *M 
O  •c 

ti  <t 

A  C  -H 
*  .£  -d 

.*’1 
e  c 

o  C 


1 


i *d 
9  9 

.•J  «• 

vi  * 

*0  .a 


O  4  A  - 
il  rt  41  C 
o 

>.  *  -  «H 

•  C  41 
>^304 
W  H  «H  44 

3  41  * 

9  0 

•  0  ,C 

4  t>  0  U 

44  0  H  « 
4  0  9 

*0  W  C 

W 
O 

.•  W 

«  *e 
£  0  P  0 
W  44  -rl 
0  0  V 

*M  «d  4J  <j 

wo  S 

&ss* 

a  *58 

n  .i  k  a 

«  .-I 

41  A  9  9 
9  *  .C  A 
X  A  A  4> 


*  J8  c 


3 

2 


W 

44 

£ 


4 

9 

H 

* 

44 

c 

o 

*H 

44 

4 

P 

C 

9 

44 

44 

4 

II 

•o 

4 

2 


W 

44 

0 

z 


o 

U3 

4 

2 

A 

44 

e 

o 


o 

44 

9 

V 

>i 

H 

c 

o 


44 


•d 

o 


44 

O 

C 

o 

83 


*-4  4 

2  0 

22 

H  4 


O 

44 


0 

43 

44 

o 

u 

c 


•  • 
C  4 

o 

44  4 

O  4 

C  A 

.  * 

3  & 

O  C 

•  o 
44  a 

4  O 

W  4 

>1  o 

4  *d 
o  p 

O  r4 

•o  o 

c 

4  -H 
•d 

O  44 

u  o 
'-*■  c 

4  4 

8  8 

*5* 

4  4 

<H  «H 
M 

•  0 

4  44 
*H  <H 

o  w 

SB  O 


c 

o 

•H  *• 

44  9 
4  W 
44  O  4 
4  <H  44 

■8  5^ 
!£  * 
w  41 

O  J  4 

*  0 


>i  w  a* 

44  44  O 
•H  O  0 

c  S.*d 

0  4  C 

*•* 

H  4  o 

0  *H  O 

woe 
44  e  o 

m. 

4  O  «4 
W  W 

0  0=0 

3  «G  *H  > 

O  44  O  44 

“*•  ‘i  § 

•  >41 

4  4  C  O 

■H  £  #  £ 
O  *H  44 

C 

e 

44  O 

C  -H  S 

0  44  4 

•H  4  0  C 
£  44  H  O 
E  4  *H  -H 
4  «H  44 
rH  4 

•  H  144 

£  4£  4 
44  44 

*  W 

4  0  4  0 

*d  4  TJ  43 

4  4  4  44 

•  O  O  O 

w  w 

_  4  o 

B  *H  8  £3 
■H  £  -H  41 
W  44  C/1 
44  44  M 

•  C  O  O 

ZHZ«H 


is 


: 

•2  4 

S3 

e2 

44 

:  i  o 

1  4  • 

I  C  *H 

I  a.o 

4  U 


I. 

4  44 

•  -H 

4 

M 

-s 


2  -o 

S 


-  4  H  4< 
•  C  W 
<H  O  O  O 
A  A  % 
4  41  S  44 
44  4  0 

owe 

O  >  4 
4  4  •  O 

C  A  >iJ3 
O  4  44 

a*  4  w 

4  c  w  c 
O  O  4  *H 
W  <H  S 
44  4 

0  4  O  C 
44  44  r-i  O 
■H  4  -r4  -H 

4  <H  44 
<H  4 

OH  041 
A  4  A  4 

41  - 43  W 

5  2  *0  * 

4  4  4  44 
0  0  0  0 

W  W 

b3  bS 

■H  £3  -H  41 
01  44  W 
44  44  W 

O  C  «  O 

%  H  *  «M 


•d 

4 

9 

W 


•H 

3 

4 
4  O 
W  fH 

4  O  «H 

M  4>W 

O  C  4 

44  II  C 

C  V  *H  o 
1/  .0  «n 

O  0*4  *4 
■H  141  4 

8- 


SI* 

b  <Mr{  M 

*  ©  _  T*  O 

S^SI* 

g§3*§ 

0*  » 


<*8 


S  8  8  * 8 

5  J3  W  V  A 

44  3  *4  4> 

O  4  O  £  4 
%  O*  4  3  O* 


i 


o  n 


„  H’S’S'g* 

h  5  ?l  41  ^ 

S  a  c  c  c 

s  -rl  M  <H  -H 
cw  h  Im  IW  Ik 


Q  Ot  U* 


44 

X3 

44 

4 

£W 


C 

O 

U 

44 

4 


8 

H 

U 

s 


H 

C4 

§ 


s  8 

(m 

s 
i 


c  e 

o  o 

■H  *rl 

4f  44 

44000000000000000000000000000000000  44 

0eceeeeeceeceecccccccedceec3ececccc  440444 

. . *s»MWMWW 

•  444444444444444444444444444444444444444444444444444444444444444444  *000  00 

SSSSSSSSSSSSSSSSSSSSSSSSSSS&SSSS&  •  j  • 

•'•••OO0d0000OO00000OOO#00000O00009  ••  >1  >1  Sw  >1 

IH'H'H’HH'H’HH'H'H^'H’H'H^I'H'H’H'H^'HH'H'H  00000 

. .  >4  w  w  w  w  w 

tC  w  W  w  M  W 

o  0  0  0  0  0 

(4 

•  *.  ”  “  “  “  “  “  ”  "  ”  “  “  “  “  "  “  “  "  "  "  "  "  “  "  "  "  "  "  §5  "  “  ”  " 

MHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH  fl  (n  H  H  H  H 

OfnhhhbihhbifNfMfHbifiMbihbilNBibilMhfiNlMlNiMlihlNhhhlMfM  W  •  b  W  h  h 

SUIfiafiBBBBBfiBBBBBBBBSflBIBBfififlBflSfiaB  9  S  S  B  0  H 

IHOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO  aooooo 

.QMWWWWWWWWWWWWWWWWWWWWWUWWWWMWWWM  2  a*  w  w  w  w 

9*  0  (x  (x  J«i  b* 

01  o  pi  «h 

29500000  00000000000000  0000000000  000  M  (ft  0  0  0  0 

(4  4444444444444444444444444444444444444444444441444444444444444444  H  4>  44  44  44 

t7*00000000000000000000000000000000  S  |X  0  0  0  0 

H  CQOQQQQQQQQQQQOQQQQQQQQQOQQQQQQQQ  C  O  Q  O  O 

*  •  o  0000  00  o  000000000  o  000  00  000000000  01  •  o  o  o  o 

ecuuoooouooooouoouuuuouuuououoouuu  teuuuu 


34 


tsm  ■  Amssw  .  >TT1  mast  txmmnt  mama  XAaa 


l~ 


A 

I 

JJ 

S  «w 

3  O 

I  s 


C  *H 

O  *0 


*H  * 

85 

•H  > 

0  • 

V  u 

0  - 
J3  « 

A  rH 
«H  • 

C  W  C 

*H  0 
tPH  «rl 
9  0  4) 
A  M  0 
P  H 

§  S 

0  H 


O 

-P 


.d 

.C  0  4) 

0* 

3  jc  d 

O  ^  rl 

"  -H  tn 

0  H 
>  o 

8  °  - 
•H  A  0 
0)  M  • 
4J  <H 
0  4)  -H 
SR  0  «M 


C  ■ 

0 

0 

0 

-G 


C 

o 

*H 

4) 

n 

H 


c 

*& 

J§ 

o 

•p 


3 

0 

%4 

5 

0 

A 


0 

0 

n 

0 

0 

3 

0 

e 


■d 

g  i .  s » 

•  «H  CL  .G  *H 

3  t>  3  41  11  3 
<M  C  K  *0  0  H 

•  4J  0  C  *H  0 

0  0  0  «3  > 

0  *0  M  p 
O  O  O  0 
O  0  hi  *  %4  A 
3  A  •  o  4) 

SSSsH  § 

•  m  m 

•H  3  • 

,13  P  -H 

0  0  0 

P  0  <C 

0  P 

0  4) 

0  O  >i 

C  A  0 

0  H 

.G  *0  *g 

4J  H  0 

3  -H  - 

P  O  <M  m  rH  A 

0  0  *H  H  * 

G  O  0  -H  ~ 

6  0  «H  » 

M  -H  £L.d  V 

O  CO  0  P  ^  CO 

P  *4  CM 

M  0  0  M  O  O  • 

0  SR  *d  o  >  * 

A  _  C  *4  p  o  at 

P  C  3  0 

0  0  O  0  C  >tO 

.G  .0  A  N  H 

»  P  -HOC  % 
0  0  A  O  0 

•  *A  A  rH 

0  =  4)  4)  0  -H 

P  C  C  P  *H  «W 

0  8  (3  0  0 

o  -H  >  jC  —  H 

■4  U  A  9  4)  <=>  O 

■9  P  •  M 

c  "0  -h  6  Sum  p 

•H  C  S  P  l)H  C 

3  8  -H  O 

*0  O  0  -H  rl  v  o 

C  jQ  -H  X  -H  O 

3  P  U  4  A  ♦  0 

O  3  0  8  o  A 
A  O  A  A  ^  43 

A  -H  0  O 

3  tM  (4  .d  P  A  C 

Oh  op  ft  0  ti 


A 

O 

M 


OIAIO0  AWH 
O  V  Ol  0t  01  o  o 

Ot«0IOIOO 

c«aa«»A 

o  o  o  o  o  o  o 


A 

o 

M 

0« 

p 

0 

c 


m  r-  0*  oo  m  o  at 

CO  A  A  A  A  O  A 
NA  AAAO  A 
o  co  oo  co  co  at  oo 

o  o  ©  o  o  o  o 


•§ 

& 

P 

0 

c 


0000000000000000000000000000 

PPPPPPPPPPPUPPPPPPMPPPPMPPMP 

^&5,&j?li’!i’&i’^5,j^SrSrS,5,S,S,S’SpS’5’5’5’S?j 

PPPPPPPPPPMPPPPPPMMPPPPPPPPP 

PMPPMPPPPPMPPPPMPPPPPPPPPPPP 

HHHHHHHHHHHHHHHHHHHHHHHHHHHH 

ftftftftftftfthlftftftftftftftklMftftlMftftftftftlwbft 

868888886888  1  B  G  fl  S  0  fl  fi  f  8  6  8  8  6  6  8 

OOOOOOOOOOOOl  coooooouooooooo 

febibjfcfolatfiibhfetelNhbbifefxfriEMfctefahjbhtetKfo 

'S'S'S'S'S'S'S'S'S-S'S'S'S'S^-S'SiS'S’S'S'S'S'S-SIg'S'S 
1 5  5  «  «  «  &  M  ?!  *  M  S:  5  &  &  3, 

4444444444444444444444444444 

aaaaaQaQaaoaooaooaaoaoooaaaa 

oooooooooooooooooooooooooooo 

uuuuuuoouuuuuauouuouuuuauouu 


o» 

o 

m 

H 

O 

at 

l 

H 

o 

I 

m 

o 

§ 


M 


CL 


M 

4) 


o  o  m  f*  in  r*  fH 

O  O  00  H  V  H  V 

o  o  o  m  vo  m  m 

OOO  ANVOO 
(0  0  MO  If)  V  t 


tH  >1  c  d  c  c  c 


CL 


o  A  A  H  H  A  A 
onaorinm 
o  at  o  m  m 
o  u>  co  at  co  ^  oo 

oo  r~  v>  ^  ^  ^ 


>i  C  d  C  C  C  C 


8* 


A 

* 


8* 


o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  ©  o  ©  o  o 
o  o  o  o  o  o  o 

o  o  o  o  o  o  o 


o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 

o  in  o  m  o  m  o 
h  n  w  f-  a 


CL 


A 

9 

rH 

-H 

(L 


CL 


M 

A 


1 

1 

H  H 

■H 

s 

0A  WMfltf  If) 
fH 

1 

H 

CD 

c 

c 

c 

o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
©  o  o  o  o  o  o 

m  m  m  m  m  m  m 

H  H  H  H  H  H  »H 


©  ©  O  ©  ©  ©  © 
O  O  O  O  O  O  O 
O  O  ©  O  ©  O  O 
O  O  O  O  O  O  O 

o  m  o  m  o  m  o 
rtnf »h a 


C 

o 


<  M  V  X 


35 


0.0000  30.0000  16  y  8.0000  0.3788 
15.0000  30.0000  5  n  7.0109  0.8991 
30.0000  30.0000  7  n  5.7309  0.9000 
45.0000  30.0000  7  n  4.6503  0.9000 


PI 

vo 

co 

• 

O 

co 

m 

o 

m 

CO 

• 

00 

O  Pi 

CM  10 

CO  00 

• 

VO 

vo  co 

CM  CM 

• 

CO 

00 

m  co 

o 

• 

O 

PI 

o 

p 

vo 

o  cn 

in 

m 

M 

O 

CO 

CM 

CM  P 

PI 

vo 

N 

p 

CM 

CM  VO 

co  in 

N 

p 

co  vn 

vo  CO 

CO  vo 

N 

CM 

PI 

vo  o 

to 

vo 

N 

O 

CM  H 

pi 

CM 

PI  vo 

pi 

in 

o 

ft 

p 

Pi 

p 

m 

CM 

CO  vo 

o  o 

p  in 

p 

CM  H 

PI  PI 

in  in 

00 

p 

p  n  vo 

m 

«rl 

o 

CO 

CO 

<p 

CM 

h  m 

o 

CO 

« 

o 

m 

vo 

CO 

pi 

Pi 

CO 

« 

CO 

co  m 

®  p 

00  00 

« 

p 

ft  u> 

pi  vo 

ao  co 

« 

P 

m 

co 

cm  p 

co 

<0 

* 

o 

Pi 

in 

•p  p 

Pl  co 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•  • 

•  • 

• 

•  • 

*  • 

•  • 

• 

• 

• 

•  • 

• 

« 

• 

• 

• 

• 

« 

•  • 

04 

CO 

vo 

CO 

CO 

CO  *0 

04 

p 

vo 

CO  lO 

CO  4* 

04 

p- 

vo  'p 

3 

3 

co  * 

n, 

p  »o  m 

p  cn 

CO 

04 

CO 

VO 

m 

' 

CO  ^ 

•  •  •  •  • 


c  c  c 

o  vo  m 


u 


>1  c  c  c  c  c  «s 
<0  o  p  p  w  o 


co  o  p  p  m  o  «h 

H  H  *  f 


e  c  c  c  c  c  c 
p  p  ft  p  vo  vo  vo 


u 

•Ai 

1 


c  c  c  c  c  c  c 

io  h  o®  r*  f  v  t* 


M 

A) 

! 


c  c  c  c  c  c  c 

IO  P  P  P  P  VO  P 


u 

■P 

1 


>t  c  c  c  c  c  c 


^  m  p*  h  p  p  o» 
H 


o 

o 

o 

c 

o 

o 

O 

o 

o 

o 

o 

C 

O 

o  o 

o 

o 

o 

o 

c 

o 

o  o 

o 

o 

o 

o 

c 

o 

o  o 

o 

o 

o 

o 

o 

0 

o 

o 

o 

o 

o 

o 

o 

0 

o 

o  o 

o 

o 

o 

o 

0 

o 

o  o 

o 

o 

o 

o 

0 

o 

O  O 

o 

o 

o 

o 

o 

fl 

o 

o 

o 

o 

o 

o 

o 

r-i 

o 

o  o 

o 

o 

o 

o 

H 

o 

o  o 

o 

o 

o 

o 

fl 

o 

o  o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

•rt 

o 

o  o 

o 

o 

o 

o 

•H 

o 

o  o 

o 

o 

o  e 

o 

o  o 

o 

o 

• 

• 

• 

a. 

• 

• 

• 

• 

• 

• 

• 

04 

• 

•  • 

• 

• 

• 

• 

04 

• 

«  • 

« 

• 

• 

• 

04 

o 

o 

o 

• 

m 

vn  m 

m  m  m  m 

• 

o 

o  o 

o 

o 

o 

o 

• 

in  m  in  m  m 

m 

in 

• 

o 

o  o 

o 

o 

tn 

Cl 

CO 

Mf 

MP 

4* 

vo 

10  vo 

vo 

VO 

vo 

vo 

p 

p  p 

p 

p 

p 

p 

Ol 

PI  A  Pl  Pl 

o  o  ©  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 

m  m  to  u)  m  m  m 
o  o  o  o  o  o  o 
H  H  H  H  »-4  H  *-4 


O  O  O 
O  O  O 
o  o  o 
o  o  o 

o  m  o 
vo  p  oi 


4J 

4 

H 

8* 


o  o  o  o  o  o  o 

o  o  o  o  o  o  o 

o  o  o  o  o  o  o 

o  o  o  o  o  o  o 

omoinomo 
r-«  cn  *p  vo  p  pi 


4> 

m 

r> I 

a, 

m 


o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 

o  m  o  in  o  in  o 


v 

* 


8* 


o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 

o  m  o  in  o  m  o 
Hwvrw 


8* 


o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 

o  m  o  in  o  in  o 
HPJvrwpoi 


n 

H 

«rl 

8* 


o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  ©  o 

o  in  o  in  o  in  o 
h  n  ^  vo  p  oi 


36 


epllat  epllan  numltec  epitlze  net  pcob 

0.0000  120.0000  17  y  8.0000  0.8222 


► 

» 


Cl 

o 

Ol 

Ol 

03 

01 

A 

ri 

o 

o 

ri 

co 

CM 

CO 

A 

4 

01 

o 

r- 

01 

CM 

ri 

A 

f H  O 

o 

O! 

01 

CM 

m 

A 

o 

4 

p“ 

o 

H 

o 

01 

Ol 

o 

o\ 

0) 

01 

O 

o 

co 

o 

O 

o 

o 

o 

0 

o 

01 

o 

o 

0V 

O 

01 

0 

10  O 

o 

01 

o  o 

01 

o 

o 

CM 

01 

01 

o 

o 

01 

oi 

o 

On 

oi 

Ol 

O 

u 

4 

o 

o 

o 

o 

o 

o 

u 

On 

01 

o 

o 

01 

O 

01 

u 

O  O 

o 

01 

o 

O 

01 

u 

o 

4 

01 

01 

o 

o 

01 

00 

Ol 

00 

00 

00 

oi 

04 

cm 

01 

01 

01 

Ol 

01 

01 

cu 

o 

00 

01 

O! 

oo 

01 

00 

0. 

o  01 

01 

00 

01 

01 

00 

cu 

o 

CM 

00 

00 

01 

01 

CO 

o 

o 

o 

o 

o 

o 

p 

o 

o 

o 

o 

o 

o 

o 

p 

o 

O 

o 

o 

o 

O 

o 

p 

o  o 

o 

o 

O 

o 

o 

p 

o 

o 

O 

o 

o 

o 

o 

• 

c 

• 

• 

c 

c 

c 

c 

o 

o 


CO 

m 

CM 

CO 

H 

CO 

• 

O 

Ol 

co  cm  m 

r* 

4 

• 

O  4 

in 

m 

00 

p*  i0 

4 

O  4 

4  m 

H 

10 

cM 

• 

O 

o 

4 

o 

o 

4 

00 

V0 

10 

CD 

CO 

00 

P* 

N 

o 

co  r* 

CM  00 

u» 

N 

O 

CM 

10 

o 

o 

4 

4 

N 

O  r-C 

4  O 

H 

H 

in 

N 

o 

o 

*0 

m 

in 

CO 

m 

10 

CO 

10 

p* 

o 

m 

P 

o 

4 

00  01 

p* 

p» 

in 

o 

CO 

o 

CM 

01 

o 

in 

O  4 

4  00 

CO 

in 

o 

o 

10 

m 

10 

10 

in 

u> 

<o 

CM 

10 

o 

00 

4 

o 

CM 

m  H 

00 

o 

00 

4 

o 

CO 

ri 

4* 

p- 

00 

4 

O  01 

4  H 

o 

H 

00 

4 

o 

o 

o 

m 

o 

H 

00 

• 

• 

«H 

• 

• 

•  • 

• 

• 

• 

P 

• 

• 

• 

• 

« 

• 

• 

*H 

•  • 

• 

• 

T< 

• 

• 

• 

• 

• 

lO 

in 

4 

CO 

4 

4 

0* 

00 

p* 

in  4 

CO 

4 

4 

& 

00 

r* 

10 

op 

CO 

4 

4 

U< 

oo  p* 

i0  m 

<0 

4 

4 

0« 

00 

00 

r- 

m 

•4 

<0 

4 

4 

• 

• 

4 

t- 


o 

o 


o 

CM 


o 

o 

o 

o 


C 

c 

C 

C 

c 

C 

U 

>i  C 

C 

c 

G 

G 

G 

H 

>t  C 

c 

c 

6 

G 

c 

U 

>i  C 

c 

c 

G 

C 

G 

H 

G 

G 

c 

G 

G 

4 

4 

4 

4 

P 

P 

P 

P 

10 

r- 

p 

p* 

r~ 

P* 

r- 

P- 

10 

01 

4 

Ol 

10 

10 

r* 

10 

H 

r- 

10 

m 

10 

01 

p* 

p* 

10 

10 

in 

H 

10 

10 

P* 

P* 

in 

1 

H 

1 

H 

H 

1 

ri 

1 

H 

ri 

H 

C 

C 

G 

C 

o 

o 

o 

o 

© 

o 

c 

O 

o 

O 

o 

o 

o 

o 

C 

O 

o 

o 

O 

o 

o 

© 

C 

O 

o 

o 

o 

o 

O 

o 

G 

O 

O 

© 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

o 

o 

0 

o 

o 

O 

O 

o 

o 

o 

0 

O 

o 

o 

o 

o 

O 

o 

O 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

ri 

o 

o 

o 

o 

o 

o 

o 

ri 

O 

o 

o 

o 

o 

o 

o 

H 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

■H 

o 

o 

o 

o 

o 

o 

o 

*rt 

o 

o 

o 

o 

©  © 

o 

■H 

o 

o 

o 

o 

o 

o 

© 

• 

» 

• 

• 

Cu 

• 

• 

• 

• 

• 

• 

• 

0< 

• 

» 

• 

• 

• 

• 

• 

o« 

• 

• 

• 

• 

• 

• 

• 

cu 

• 

* 

• 

o 

o 

o 

o 

o 

o 

4 

m 

m 

in 

m 

m 

in 

m 

4 

o 

o 

o 

o 

o  o 

o 

4 

m  m  m  m 

in  m 

m 

4 

o 

o 

o 

o 

o 

o 

o 

CM 

CM 

CM 

CM 

CM 

CM 

co 

<0 

co 

co 

<0 

<0 

co 

m 

m 

w  in 

in 

m 

in 

10 

10 

10  10 

10 

10 

10 

00 

00 

00 

00 

00 

00 

00 

H 

H 

H 

H 

H 

H 

H 

H 

r*4 

H 

H 

ri 

H 

ri  H 

ri 

ri 

ri 

ri 

H 

H 

H 

H 

H 

ri 

H 

ri 

H 

H 

rM 

H 

H 

H 

H 

o  o  o  o  o  o 
o  o  o  o  o  o 
o  o  o  o  o  o 
o  o  o  o  o  o 

momomo 
h  co  4*  10  p-  o» 


P 

4 

H 

H 

• 


o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 

o  in  ©  m  o'  m  o 


* 


o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 

o  in  o  in  o  in  o 
h  co  10  r-  01 


p 

4 

H 

■H 

A* 


o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  o  o  o 
o  o  o  o  ©  o  o 

omdmoioo 

h  <o  4  10  p* 


p  oooooooz 

4  00000004$ 

H  OOOOOOOZ 

■H  OOOOOOO 

P<  . P 

•  omomomoo 

w  <o  4  10  p*  01  * 


S 


37 


OUTPUT  DATA 


NetSim  writes  two  output  data  files;  the  graphics  file  (.gra)  and  the  list  file  (.1st).  The 
formats  of  these  files  are  described  in  detail  in  Appendix  A.  The  graphics  file  from 
the  sample  NetSim  run  is  given  in  Table  4.6.  This  includes  the  data  in  the  control  file 
and  the  data  needed  to  make  contour  plots  (station  names  and  coordinates,  detection 
thresholds,  probabilities,  and  location  uncertainties).  The  last  nine  columns  of  the 
matrix  in  Table  4.6  are  zero  because  these  are  location  uncertainties  and  the  sample 
NetSim  run  only  calculated  detection  thresholds.  Table  4.7  lists  the  first  three  pages  of 
the  list  file,  and  one  page  for  epicenter  33  (there  is  one  page  of  detection  parameters 
for  each  ep;center). 


SNAPCON  TRANSCRIPTION 

A  complete  transcription  of  an  interactive  SnapCon  session  to  generate  a  contour  map 
of  the  detection  thresholds  calculated  in  the  sample  NetSim  run  starts  on  page  46  of 
this  report.  Comments  are  separated  from  the  transcription  by  a  vertical  line  on  each 
page.  This  transcription  is  included  as  an  example  of  the  SnapCon  user-interface. 
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Table  4.6.  Sample  output  graphics  file  (continued). 
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Table  4.7.  Sample  output  list  file. 
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Table  4.7.  Sample  output  list  f.te  (continued). 
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Table  4.7.  Sample  output  iist  file  (continued). 
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user  selects  "d"  to  modify  contour  options. 
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DETECTION  THRESHOLDS 


Figure  4.6a  plots  contours  of  the  detection  thresholds  calculated  for  the  NetSim  sample 
run  at  1.0  Hz  (this  is  the  result  of  the  SnapCon  session  described  above).  The  detec¬ 
tion  threshold  is  MLg  3.7-4.2  for  epicenters  in  the  Soviet  Union.  Figure  4.6 b  plots  the 
detection  thresholds  a t  10  Hz.  These  results  were  obtained  by  changing  the  frequency 
in  the  control  file  (Table  4.1)  to  10  Hz,  and  repeating  the  NetSim  calculations 
described  for  the  first  sample  run.  The  90%  MLg  threshold  for  detecting  3  P  waves  at 
10  Hz  is  2. 8-5.3  for  epicenters  in  the  Soviet  Union.  The  threshold  at  10  Hz  is  lower 
than  the  threshold  at  1  Hz  for  epicenters  in  areas  with  dense  station  coverage,  but  the 
threshold  is  much  higher  at  10  Hz  than  at  1  Hz  for  epicenters  in  areas  with  sparse  sta¬ 
tion  coverage.  This  is  because  the  frequency  of  the  maximum  signal-to-noise  ratio 
(SNR)  decreases  with  increasing  epicentral  distance.  It  is  clear  from  these  examples 
that  frequency  dependent  input  parameters  are  required  to  obtain  accurate  estimates  of 
detection  capability  at  regional  distances. 

4.2  Frequency-Dependent  Input 

This  section  presents  the  results  of  a  sample  NetSim  run  that  includes  frequency  depen¬ 
dence  in  the  source,  attenuation,  station  noise,  and  beam  gain.  We  use  the  same  input 
data  files  that  were  described  in  Section  4.1,  except  that  we  have  modified  the  control 
file  to  include  11  frequencies  between  1.0  and  20.0  Hz  (Table  4.8).  The  results  of  the 
calculation  are  plotted  in  Figure  4.7.  The  frequency-dependent  90%  MLg  threshold  for 
detecting  3  P  waves  is  2.6-3.4  for  epicenters  in  the  Soviet  Union.  These  thresholds 
are  less  than  the  thresholds  for  either  1  Hz  or  10  Hz  (Figure  4.6).  This  is  because  the 
frequency  of  the  maximum  SNR  is  greater  than  1  Hz  and  less  than  10  Hz  for  epicen¬ 
tral  distances  important  for  determining  the  detection  threshold  for  the  network.  For 
example.  Figure  4.8  plots  the  frequency  of  the  maximum  SNR  as  function  of  epicentral 
distance  for  an  epicenter  in  western  USSR  (60°  N,  60°  E).  Table  4.9  lists  the  detec¬ 
tion  parameters  for  this  epicenter.  The  frequency  of  the  maximum  SNR  is  4-7  Hz  for 
all  stations  that  have  a  probability  of  detection  greater  than  zero. 

4.3  Probability-Weighted  and  Monte  Carlo  Location 

There  are  two  ways  that  NetSim  can  account  for  the  effect  of  undetected  phases  on 
location  uncertainties;  the  probability-weighted  approach  (used  in  SNAP/D)  and  the 
Monte  Carlo  approach  described  in  Section  2.3.2  In  this  section  we  compare  the  loca¬ 
tion  uncertainties  computed  using  these  two  methods  for  events  near  the  3  P  detection 
threshold  of  the  network.  We  also  compare  the  two  approaches  for  one  epicenter  in 
western  USSR  as  a  function  of  event  size. 

We  use  the  same  input  data  files  for  our  location  uncertainty  simulations  that  we  used 
for  our  detection  threshold  simulations  (Section  4.3).  The  only  difference  is  the  runtyp 
parameter  in  the  control  file  in  Table  4.8  was  changed  to  detloc  to  compute  the  median 
and  90th  percentile  Monte  Carlo  location  uncertainties.  The  probability-weighted 
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CONTROL  DATA 

CONTROL  riLB :  . / fdrunl • out 

Frequency-Dependent  Detection  Threihold* 

PDrunl 

./results 

./ 

detection 

threshold 

<P/3) 

P 

alg 

3.000  1.000  8.000  0.900 

0.900  100  0.500  0.900 

fixed  0.000 
11 

1-  2.  3.  4.  5. 

6.  7.  8.  10.  15.  20. 
eurss.epi 
granite .  end 
Tranlte . sgr 
able. pod 
'ble.pgr 
juOxi?.  sta 


0.900  20.000 


Table  4.8.  Sample  control  file  (frequency-dependent  input). 


51 


XS  mZZ of  ,he  90,6 


vJV/iw,UMg 


Frequency 


Frequency  of  Maximum  SNR  (P  Wave) 


Epicentral  Distance  (Degrees) 


Figure  4.8.  Frequency  of  the  maximum  SNR  as  a  function  of  epicentral  distance  for  P 
phases.  This  figure  plots  the  values  in  Table  4.9  for  an  epicenter  in  the  western  USSR 
(60°  N,  60°  E). 
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output  List  FILE:  ./r*«ult./FDrunl.l»t 


EPICENTER 


STATION 

1  NOR 

2  SBIO 

3  GUPO 

4  GOB 

5  RSTN 

6  ANTO 

7  MAN 

8  MAJO 

9  BRW 

10  TOL 

11  PKR 

12  FCC 

13  NR2 

14  SI 

15  S2 

16  S3 

17  S5 

18  S6 

19  S7 

20  S8 

21  S9 

22  S10 

23  Sll 

24  S12 

25  S13 

26  S14 

27  SIS 

28  S16 

29  SI 7 

30  S18 

31  S19 

32  S20 

33  S21 


33 


1AT-  60.000  LONG-  60.000 

EVENT  SIZE 

NUMBER  OF  ITERATIONS 

SEMI -MAJOR  AXIS  (LO) 

SEMI -MINOR  AXIS  (LO) 

STRIKE  (LO) 

DEPTH-INTERVAL  (LO)  : 

MONTE  CARLO  FRACTION  : 


2.598  NETWORK  PROBABILITY  (%) 

6  ITERATION  BOUNDS  EXCEEDED 

0.000  SEMI -MAJOR  AXIS  (HI) 

0.000  SEMI -MINOR  AXIS  (HI) 

0.000  STRIKE  (HI) 

0.000  DEPTH-INTERVAL  (HI) 

0.000 


89.977 

n 

0.000 

0.000 

0.000 

0.000 


ANC. 

DIST. 


P  -NAVE 

STATION  FREQ.  OF  AHPL. 
AZM  PROB.  MAX.  S/N  ATTBN.SD 


23.424 

293.015 

40.822 

133.209 

30.176 

167.990 

42.343 

331.168 

57.450 

357.039 

26.282 

232.411 

63.054 

108.215 

53.293 

78.814 

46.220 

15.405 

43.681 

272.471 

56.361 

16.152 

59.520 

344.754 

16.995 

321.538 

30.040 

86.255 

49.182 

26.213 

45.365 

63.710 

25.076 

185.806 

15.000 

180.273 

17.406 

142.862 

17.500 

217.777 

7.765 

234.812 

6.515 

136.684 

15.069 

97.557 

24.486 

81.393 

33.605 

69.290 

48.806 

48.811 

5.503 

337.414 

11.081 

54.042 

22.221 

56.415 

32.829 

48.787 

42.002 

38.251 

16.675 

5.690 

22.797 

23.309 

0.000 

2.000 

0.000 

1.000 

0.000 

2.000 

0.000 

1.000 

0.000 

1.000 

0.000 

2.  COO 

0.000 

1.000 

0.000 

1.000 

0.000 

1.000 

0.000 

1.000 

C.000 

1.000 

0.000 

1.000 

0.000 

3.000 

0.000 

2.000 

0.000 

1.000 

0.000 

1.000 

0.000 

2.000 

0.002 

4.000 

0.000 

3.000 

0.000 

3.000 

0.888 

5.000 

0.984 

7.000 

0.001 

4.000 

0.000 

2.000 

0.000 

2.000 

0.000 

1.000 

0.999 

7.00C 

0.210 

4.000 

0.000 

2.000 

0.000 

2.000 

0.000 

1.000 

0.000 

4.000 

0.000 

2.000 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 

0.300 


Table  4.9.  Sample  list  file  (frequency-dependent  input). 
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location  uncertainties  were  calculated  by  changing  nwncarlo  in  the  control  file  to  1 
(see  Appendix  A).  Figure  4.9  plots  the  probability-weighted  estimate  of  the  semi¬ 
major  axis  (SMA)  of  the  epicentral  error  ellipse.  The  probability-weighted  SMA  is 
between  14  and  40  km  for  epicenters  in  the  Soviet  Union.  The  results  for  the  median 
and  90th  percentile  Monte  Carlo  realizations  are  plotted  in  Figure  4.10.  The  median 
Monte  Carlo  SMA  is  between  13  and  45  km  for  events  in  the  Soviet  Union.  This  is 
similar  to  the  range  produced  by  the  probability-weighted  approach  (the  average 
difference  between  the  median  Monte  Carlo  SMA  and  the  probability-weighted  SMA 
is  3.3  km  for  epicenters  in  the  Soviet  Union).  The  90th  percentile  Monte  Carlo  reali¬ 
zation  SMA  is  between  20  and  235  km  for  events  in  the  Soviet  Union,  which  is 
significantly  larger  than  the  median  Monte  Carlo  SMA  or  the  probability-weighted 
SMA. 

The  results  in  Figures  4.9  and  4.10  are  for  events  at  the  90%  MLg  threshold  for  detect¬ 
ing  3  P  phases.  Since  azimuth  and  arrival  time  can  be  estimated  for  each  P  detection, 
there  should  be  more  than  enough  data  to  constrain  the  event  epicenter  (assuming  fixed 
depth).  We  found  that  at  least  95%  of  the  Monte  Carlo  realizations  had  enough  data 
to  constrain  the  event  epicenters  in  the  Soviet  Union.  However,  we  also  found  that  the 
probability-weighted  approach  increasingly  underestimates  the  median  Monte  Carlo 
SMA  as  the  event  magnitude  decreases.  For  example,  Figure  4.11  plots  histograms  of 
the  SMA  computed  using  the  Monte  Carlo  approach  for  three  different  event  magni¬ 
tudes  at  an  epicenter  in  western  USSR  (60°  N,  60°  E).  One  hundred  realizations  were 
generated  for  each  magnitude.  Figure  4.11a  is  for  MLg  3.0,  which  is  greater  than  the  3 
P  detection  threshold  at  this  epicenter.  The  probability-weighted  SMA  is  nearly  equal 
to  the  median  Monte  Carlo  SMA  (19  km),  and  less  than  the  90th  percentile  Monte 
Carlo  SMA  (27  km).  Figure  4.11b  is  for  MLg  2.58,  which  is  equal  to  the  3  P  detec¬ 
tion  threshold.  Again  the  probability-weighted  SMA  (25  km)  is  close  to  the  median 
Monte  Carlo  SMA  (27  km),  but  it  is  substantially  less  than  the  90th  percentile  Monte 
Carlo  SMA  (140  km).  Figure  4.11c  is  for  MLg  2.23,  which  is  equal  to  the  2  P  detec¬ 
tion  threshold.  In  this  case,  only  79%  of  the  Monte  Carlo  realizations  had  enough  data 
to  constrain  the  location.  The  probability-weighted  SMA  is  34.6  km,  while  the  median 
and  90th  percentile  Monte  Carlo  estimates  are  both  140  km.  It  is  clear  from  Figure 
4.11  that  a  decrease  in  the  event  magnitude  results  in  an  increase  in  the  discrepancy 
between  the  probability-weighted  and  Monte  Carlo  approaches. 

The  probability-weighted  approach  is  most  accurate  when  an  individual  datum  is  not 
very  important  for  constraining  the  location  [Ciervo  et  al. ,  1985].  Obviously,  when 
few  data  are  available  to  locate  an  event  (e.g.,  for  most  events  near  the  detection  thres¬ 
hold),  each  datum  is  very  important  It  is  likely  that  the  probability-weighted  approach 
underestimates  location  uncertainties  for  small  events  because  data  with  low  probabil¬ 
ity  of  detection  provide  unique  information  that  is  important  to  constraining  location. 
It  appears  that  the  probability-weighted  approach  does  not  adequately  reduce  the  con¬ 
tribution  of  these  data  to  the  location  uncertainty. 
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RUN  OflTE 

C0NT0UR 


3/  1/1990  fiT  17.  7.38 

FR0M  10.000  T0  110.00 


C0NT0UR  INTERVAL  0F  20.000 


Figure  4.9.  Contours  of  the  semi-major  axis  (SMA,  in  kilometers)  computed  using  the 
probability-v-eighted  approach  for  events  at  the  3  P  detection  threshold.  Uncertainties 
are  given  at  the  90%  confidence  level.  Depth  is  constrained  to  the  surface. 
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SEMI-MRJ0R  UNCERTfllNTr  AXIS  IKHBNE  FILE  t*  SEMI -MRJ0R  UNCERTflINTr  AXIS  1KMBNE  FILE 


Semi-Major  Axis  (km) 

Figure  4.11.  Histogram  of  location  uncertainties  (SMA)  determined  from  100  Monte 
Carlo  realizations  of  location  solutions  in  the  western  Soviet  Union  (60°  N,  60°  E). 

The  probability-weighted  (PW),  the  median  Monte  Carlo  realization  (MC50),  and  the 

90th  percentile  Monte  Carlo  realization  (MC90)  SMA  values  are  marked.  The  fraction 

of  Monte  Carlo  realizations  with  enough  data  to  constrain  the  epicenter  is  also  listed.  t 

Event  magnitudes  are  (a)  MLg  3.0,  ( b )  MLg  2.58  (3  P  detection  threshold),  and  (c)  MLg 

2.23  (2  P  detection  threshold). 
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5.  SUMMARY 


This  report  describes  a  computer  program  called  NetSim  that  simulates  the  detection 
and  location  capability  of  seismic  networks  that  include  regional  stations  and  arrays. 
The  main  advantage  of  NetSim  over  other  simulations  programs  is  the  incorporation  of 
frequency  dependence  in  the  estimates  for  the  source,  attenuation,  station  noise,  and 
array  gain.  This  is  very  important  for  regional  networks,  as  is  demonstrated  by  the 
comparison  of  frequency-dependent  and  frequency-independent  detection  thresholds  for 
a  network  of  33  stations  in  and  around  the  Soviet  Union.  The  results  are  shown  in 
Table  5.1. 


Table  5.1  The  90%  MLg  thresholds  for  a  network  of  33  stations. 

center  expand  tab  (:)  ;  c  c  c  c.  _  FigureiDetection:  :90%  MLg 

Numher:Criterion:Frequency:Threshold  _ 

4.5a: 3  PA  Hr.  (fi\ed):3.7-4.2  4.56:3  />:10  Hz  (Sxed):2. 8-5.3  4.6:3  PA-20  Hz  (vary- 
ing):2.6-3.4 


These  simulations  are  based  on  regional  wave  attenuation  and  noise  conditions  at 
NORESS,  where  the  frequency  of  the  maximum  SNR  decreases  with  increasing  epi- 
central  distance.  The  frequency-dependent  detection  thresholds  are  lower  than  the 
thresholds  at  either  1  Hz  or  10  Hz  because  the  frequency  of  the  maximum  SNR  at 
ranges  that  control  the  detection  threshold  is  typically  between  3  and  7  Hz.  It  is  clear 
from  this  example  that  the  detection  capability  of  regional  networks  cannot  be  simu¬ 
lated  accurately  unless  frequency-dependent  factors  are  included. 

Other  important  advantages  of  NetSim  over  similar  computer  programs  are:  (1)  noise 
levels  for  secondary  phases  are  based  on  the  coda  of  previous  arrivals  (and  depend  on 
event  magnitude  and  epiceutral  distance),  (2)  arrival  time  and  azimuth  standard  devia¬ 
tions  are  used  to  estimate  location  uncertainties,  and  (3)  the  effect  of  undetected  phases 
on  location  uncertainties  can  be  estimated  using  a  Monte  Carlo  approach  that  mimics 
the  way  location  uncertainties  vary  in  practice.  The  probability-weighted  approach 
(used  in  SNAP/D)  assumes  that  all  data  contribute  to  the  location,  but  the  standard 
deviations  are  divided  by  the  square  root  of  the  probability  of  detection.  Thus,  arrival 
times  and  azimuths  £om  phases  with  low  probability  of  detection  will  not  contribute 
much  to  the  location  solution  since  their  effective  standard  deviations  are  large.  How¬ 
ever,  comparison  to  results  from  a  more  realistic  Monte  Carlo  approach  suggests  that 
the  probability-weighted  uncertainties  underestimate  the  true  location  uncertainties  for 
small  events.  This  implies  that  the  probability-weighted  approach  doe  .  not  adequately 
degrade  the  contribution  of  data  with  low  probability  of  detection. 
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APPENDIX  A:  INPUT/OUTPUT  DATA 


This  appendix  describes  the  NetSim  input  and  output  data  and  file  formats.  The  input 
to  NetSim  is  read  from  standard  ASCII  files.  These  input  files  can  be  updated  interac¬ 
tively,  and  new  files  can  be  written  using  the  NetSim  data  handling  module.  NetSim 
output  is  also  written  to  standard  ASCII  files.  One  of  these  output  files  is  read  by 
SnapCon,  which  plots  the  results  (detection  thresholds,  probabilities  for  fixed  event 
size,  or  location  uncertainties)  as  contours  on  regional  or  world  maps. 

There  are  five  categories  of  input  data:  control,  source,  propagation,  site/station,  and 
noise  (see  Section  3).  Each  category  has  one  or  more  types  of  input  data  files.  We 
identify  these  different  data  files  by  a  three-letter  suffix  appended  to  the  data  file  name 
(these  suffixes  are  included  in  tables  in  this  appendix  for  each  input  data  file).  The 
NetSim  computer  program  does  not  impose  this  file  naming  convention,  but  we  have 
found  it  to  be  convenient.  Tables  A.1-A.15  list  the  format  of  the  input  and  output 
data  files.  Indented  lines  indicate  a  continuation  from  the  previous  line.  Sample  files 
are  given  in  Section  4. 

NetSim  includes  provisions  for  up  to  ten  seismic  waves  (or  phases).  The  wave  input  is 
ordered  such  that: 


Wave(l) 

— 

P  (or  Pn) 

Wave(2) 

- 

pP 

Wave(3) 

= 

Pg 

V/ave(4) 

= 

S  (or  Sn) 

Wave(5) 

= 

Lg 

Wave(6) 

= 

LR 

Wave(7) 

= 

01 

Wave(8) 

= 

02 

Wave(9) 

= 

03 

Wave(10) 

= 

04 

where  the  last  four  are  slots  for  additional  seismic  phases. 


INPUT  DATA 

Control  Data.  The  control  data  define  the  nature  of  the  NetSim  run  (detection  or  loca¬ 
tion  or  both,  network  detection  probability,  detection  criteria,  etc).  There  is  only  one 
input  file  in  this  category;  the  control  file.  Table  A.l  shows  the  format  of  the  control 
file.  The  following  list  defines  each  input  parameter. 
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Table  A.l.  Control  File 


CONTROL  FILE 

Suffix:  .cnt 


title 

outpre 

outdir 

dir 

runtype 

subtype 

detcrit 

locwavs 

sizetype 

evsize  evmin  evmax  probdet  cutlo  regdist 
probloc  numcarlo  prctlo  prcthi 
depthfix  eventdepth 
numrunfrq 

runfrq(l)  ...  runfrqf numrunfrq ) 

file.epi 

file.smd 

file.sgr 

file.pmd 

file.pgr 

file.sta 
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Control  File  (.cnt) 


1.  title  -  Title  of  the  NetSim  run. 

2.  outpre  -  Prefix  for  NetSim  output  files.  The  output  files  are  named  ontpre. gra 
and  outpre.lst  (these  are  described  below). 

3.  outdir  -  Directory  for  NetSim  output  files. 

4.  dir  -  Directory  for  the  NetSim  input  files. 

5.  runtype  -  NetSim  nm  type.  Acceptable  values  are  "detection,"  "location,"  or 
"detloc"  (detection  and  location). 

6.  subtyp  -  NetSim  run  subtype.  Acceptable  values  are  "threshold"  and  "probabil¬ 
ity." 

7.  detcrit  -  Detection  criteria.  The  format  is  the  same  as  that  described  in  the 
SNAP/D  user’s  manual  ( Ciervo  et  al.,  1985).  Two  examples  are: 

a.  P/3;  P-wave  detection  by  at  least  three  stations. 

b.  (P/3  +  (P/2  *  Lg/1)  +  (P/1  *  Lg/2)  +  Lg/3);  at  least  three  P  or  Lg  waves 
detected  by  the  network.  Literal  translation:  P-wave  detection  at  three  sta¬ 
tions,  or  P- wave  detection  at  two  stations  and  Lg- wave  detection  at  one  sta¬ 
tion,  or  P- wave  detection  at  one  station  and  Lg- wave  detection  at  two  sta¬ 
tions,  or  Lg- wave  detection  at  three  stations. 

8.  loewavs  -  Waves  to  use  in  location  simulation  (these  may  be  different  from  the 
waves  used  for  detection). 

9.  sizetype  -  Event  size  type.  Acceptable  values  are  "moment,”  "yield,"  "mb,"  "Ms," 
or  "mLg."  If  "moment"  is  specified,  then  evsize,  evmin,  and  evmax  are  specified  as 
log  moment  (nt-m).  If  "yield"  is  specified,  then  evsize,  evmin,  and  evmax  are 
specified  as  log  yield  (kilotons). 

10.  evsize  -  Event  size  to  use  in  a  probability  run  or  in  a  location  run  with  fixed 
event  size.  The  units  of  evsize  are  determined  by  sizetype. 

11.  evmin  and  evmax  -  Minimum  and  maximum  event  size  to  use  in  a  detection 
threshold  run.  The  units  are  determined  by  sizetype. 

12.  probdet  -  Probability  of  detection  for  the  network. 

13.  cutlo  -  Low-probability  cut  off.  If  detection  probability  is  less  than  cutlo  for  any 
wave  at  any  station,  then  this  probability  is  set  to  zero. 

14.  regdist  -  Division  between  regional  and  teleseismic  distances  (in  degrees). 

15.  probloc  -  Confidence  level  for  location  uncertainty. 

16.  nwncarlo  -  Number  of  Monte  Carlo  location  iterations  (=1  for  probability- 
weighted  location  uncertainties). 

17.  prctlo  -  Percentage  of  Monte  Carlo  location  runs  with  uncertainty  less  than  the 
lowest  of  the  two  uncertainties  reported  in  the  output  file  (e.g.,  the  median  Monte 
Carlo  realization  has  prctlo  =  0.50). 

18.  prethi-  Percentage  of  Monte  Carlo  location  runs  with  uncertainty  less  than  the 
largest  of  the  two  uncertainties  reported  in  the  output  file  (e.g.,  the  90th  percentile 
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Monte  Carlo  realization  has  prcthi  =  0.90). 

19.  depthfix  -  Event  depth  constraint  for  location  uncertainties.  Acceptable  values  are 
"fixed"  and  "free." 

20.  eventdepth  —Event  depth  in  kilometers. 

21.  numrunfreq  -  Number  of  frequencies  to  use  in  the  simulation. 

22.  runfrq(i),  i  =  1,  numrunfreq  -  Frequencies  to  use  in  the  simulation  (Hz). 

23.  file.epi  -  File  name  for  epicenter  input. 

24.  file.smd  -  File  name  for  source  media  input. 

25.  file.sgr  -  File  name  for  source  grid  input. 

26.  file.pmd  -  File  name  for  path  media  input. 

27.  file. pgr  -  File  name  for  path  grid  input. 

28.  file.sta  -  File  name  for  station  data  input 

Source  Data.  The  source  data  define  source  spectral  amplitude  and  scaling  relations 
as  a  function  of  source  medium  (e.g.,  granite,  salt,  tuff).  There  are  four  input  file 
types  in  this  category:  epicenter  grid  file,  source  medium  file,  source  medium  grid  file, 
and  the  source  spectra  files.  The  first  three  of  these  are  indexed  in  the  control  file. 
The  source  spectra  files  are  indexed  in  the  source  medium  file.  There  is  a  separate 
source  spectra  file  for  each  source  medium.  Tables  A.2-A.5  show  the  format  of  the 
four  source  input  files.  The  following  list  defines  each  input  parameter. 

Epicenter  Grid  File  (.epi) 

1.  title  -  Title  of  the  epicenter  grid  file. 

2.  numepilat  numepilon  -  Number  of  latitudes  and  number  of  longitudes  in  the  epi¬ 
center  grid. 

3.  epilatl  epilonl  -  Initial  latitude  and  initial  longitude. 

4.  epilatdel  epilondel  -  Latitude  increment  and  longitude  increment. 

Source  Medium  File  (.smd) 

1.  title  -  Title  of  the  source  medium  file. 

2.  dir  -  Directory  for  source  spectra  files  (.sor). 

3.  numsormed  -  Number  of  different  source  media. 

4.  sormed(i),  i  =  1,  numsormed  -  Source  media  names  (e.g.,  granite,  tuff,  salt). 

5.  tbl.sor(i),  i  =  1,  numsormed  -  File  name  of  the  source  spectra  file. 

6.  smdrho(i),  i  =  1,  numsormed  -  Density  of  the  source  medium  (kg/m3). 

7.  ayield(i)  and  byield(i),  i  =  1,  numsormed  *•  Slope  and  intercept  of  log  moment 
versus  log  yield:  log  moment  =  ayield  x  log  yield  +  byield.  Moment  is  in  nt-m 
and  yield  is  in  kilotons. 
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Table  A.2.  Epicenter  File 


EPICENTER  FILE 

Suffix:  .epi 

title 

numepilat  numepilon 

epilatl  epilonl 

epilatdel  epilondel 

69 


Table  A.3.  Source  Media  File 


SOURCE  MEDIA  FILE 

Suffix:  .smd 


title 

dir 

numsormed 

sormed(l) 

tbl.sor(l) 

smdrho(l) 

ayield(l)  byield(i)  amb(l)  bmb(l) 
ams(l)  bms(l)  amlg(l)  bmlg(l) 
sorwavfac(l.l) ...  sorwavfac(lO.l) 
smdvel(l.l)  ...  smdvel(10,l) 


sormed(  numsormed ) 
tbl.sorf  numsormed ) 
smdrhof  numsormed ) 

ayield( numsormed)  byield( numsormed)  amb( numsormed)  bmb( numsormed) 
ams( numsormed)  bms( numsormed)  amlg( numsormed)  bmlg(numsormed) 
sorwavfac(l .numsormed) ...  sorwayfac(10, numsormed) 
smdvel(l, numsormed) ...  smdvel(  10, numsormed) 


Table  A.4.  Source  Grid  File 


SOURCE  GRID  FILE 

Suffix:  .sar 


title 

numsgrlat  numsgrlon 
sgrlatl  sgrlonl 
sgrlatdel  sgrlondel 
sormednamef  1 )  sgrlatf  1 )  sgrlonf  1 ) 


sormednamef numsgrlat  x  numsgrlon)  sgrlat(numsgrlat)  sgrlonf numsgrlon) 


Table  A.5.  Source  Spectra  File 


SOURCE  SPECTRA  FILE 

Suffix:  .sor 

title 

numsorfrq 

sorfrq(l) ...  sorfrq(numsorfrq) 
numsormom 

sormom(l) ...  sormomf numsormom) 

Source  amplitudes  for  frequency  =  sorfrq(l)  <=  comment  line 
soramp(l,l) ...  sorampf numsormom,  1 ) 

Source  amplitudes  for  frequency  =  sorfrq(numsorfrq)  <=  comment  line 
sorampfl, numsorfrq) ...  sorampf numsormom, numsorfrq) 


\ 

i 


8.  amb(i)  and  bmb(i),  i  ~  1,  numsormed  -  Slope  and  intercept  of  log  moment  versus 
mb:  log  moment  =  ami'  x  mb  +  bmb.  Moment  is  in  nt-m. 

9.  ams(i)  and  bms(i),  i  =  1,  numsormed  -  Slope  and  intercept  of  log  moment  versus 
Ms:  log  moment  =  ams  xM,  +  bms.  Moment  is  in  nt-m. 

10.  amlg(i)  and  bmlg(i),  i  =  1,  numsormea  -  Slope  and  intercept  of  log  moment 
versus  mtg:  log  moment  =  amlg  x  mlg  +  bmlg.  Moment  is  in  nt-m. 

11.  sorwayfac(j,i),  j  -  1,  10  and  i  =  1  numsormed  -  Source  excitation  factor  (k*  in 
Section  2.1). 

12.  smdvel(j,i),  j  =  1,  10  and  i  =  1,  numsormed  -  Source  m  Hum  velocity  (m/s). 

Source  Grid  File  (.sgr) 

1.  title  -  Title  of  the  source  grid  file. 

2.  numsgrlat  numsgrlon  -  Number  of  latitudes  and  number  of  longitudes  in  the 
source  medium  grid. 

3.  sgrlatl  sgrlonl  -  Southernmost  latitude  and  westernmost  longitude  in  the  source 
medium  grid. 

4.  sgrlatdel  sgrlondel  -  Latitude  increment  and  longitude  increment  for  the  source 
medium  grid. 

5.  sormedname(l),  1=1,  numsgrlat  x  numsgrlon  -  Source  medium  name  (e.g.,  gran¬ 
ite,  tuff,  salt)  for  each  grid  clement. 

6.  sgrlat(i)  and  sgrlon(j),  i  =  1,  numsgrlat  and  j  =  1,  numsgrlon  -  Latitude  and 
longitude  of  the  southwest  comer  of  each  grid  element. 

Source  Spectra  File  (.sor) 

1.  title  -  Title  of  the  source  spectra  file. 

2.  numsorfreq  -  Number  of  frequencies. 

3.  sorfrq(i),  i  ~  1,  numsorfreq  -  Frequencies  (Hz). 

4.  numsormom  ■■  Number  of  scalar  seismic  moments. 

5.  sormom(j),  j  -  1,  numsormom  -  Log  seismic  moment  (nt-m). 

6.  soramp(j,i),  j  =  1,  numsormom  and  i  =  1,  numsorfreq  -  Log  source  spectral 
amplitude  (nt-m). 

Propagation  Data.  The  propagation  data  define  tne  attenuation  and  travel  time  as  a 
function  of  path  medium  (e.g.,  stable  or  tectonic).  There  are  four  input  file  types  in 
this  category:  path  medium  file,  path  medium  grid  file,  attenuation  files,  and  travel 
time  files.  The  first  two  of  these  are  indexed  in  the  control  file.  The  attenuation  and 
trav.*.i  L.ne  files  are  indexed  in  the  path  medium  file.  There  are  separate  attenuation 
and  travel  time  files  for  each  path  medium.  Tables  A.6-A.9  s’’ow  the  format  of  the 
four  propagation  input  files.  The  following  list  defines  each  input  parameter. 
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Table  A.6.  Path  Media  File 


PATH  MEDIA  FILE 

Suffix:  .pmd 

title 

dir 

tlen(l)  gvell(l)  gvel2(l) 


tlen(10)  gvell(lO)  gvel2(10) 

numpthmed 

pthmedfj , 

tbl.atn(  1,1)  tbi.timf  1,1) 
atnstd(l.l)  sorcor(l,l)  reccor(l,l) 
timsdf  1,1)  azsd(  1,1) 


tbl.atn(10,l)  tbl.tim(10,l) 
atnsid(lOJ)  sorcor(lOJ)  reccor(lO.l) 
timsd(10,l)  azsd(10,l) 


pthmed(  numpthmed ) 

tbl.atn(l , numpthmed)  tbl.tim(  1  .numpthmed) 
atnstd(l .numpthmed)  sorcorf  1 , numpthmed)  reccor(l, numpthmed) 
timsd(l .numpthmed)  azsd(l, numpthmed) 


tbl.atn(  10, numpthmed)  tbl.tim(  10, numpthmed) 
atnstd(  10,  numpthmed)  sorcor(10, numpthmed)  reccorf  10, numpthmed) 
timsd( 10, numpthmed)  azsd(  10, numpthmed) 


Table  A.7.  Path  Grid  File 


PATH  GRID  FILE 

SuQix:  .pgr 

title 

numpgrlat  numpgrion 
pgrlatl  pgrlonl 
pgrlatdJ  pgrlondel 
pthmedname(l)  pgrlat(l)  pgrlon(l) 


pthmedname(numpgrlat  x  numpgrion)  pgrlat(numpgrlat)  pgrlon( numpgrion) 
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Tabic  A.8.  Attenuation  File 


ATTENUATION  FILE 

Suffix:  .atn 


title 

numatnfrq 

atnfrq(l) ...  atnfrqf numatnfrq) 
numatndst 

atndst(l) ...  atndst( numatndst) 

attenuation  for  frequency  =  atnfrq(l)  <=  comment  line 
atnamp(l,l) ...  atnamp(numatndst,l) 


attenuation  for  frequency  =  atnfrqf numatnfrq)  «=  comment  line 
atnamp(l, numatnfrq) ...  atnamp( numatndst, numatnfrq) 


Table  A.9.  Travel  Time  File 


TRAVEL  TIME  FILE 

Suffix:  .tint 


title 

numtimdph 

timdph(l)  ...  timdpk( numtimdph ) 
numtimdst 

timdst(l)  ...  timdst(numtimdst) 

travel- time  for  depth  =  depth(l )  <t=  comment  line 

tratim(l,l) ...  tratim( numtimdst, 1 ) 


travel-time  for  depth  =  depthf numtimdph)  <=  comment  line 
tratim(l, numtimdph) ...  tratim( numtimdst, numtimdph) 
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Path  Medium  File  (.pmd) 


1.  title  ~  Title  of  the  path  medium  file. 

2.  dir  -  Directory  for  the  attenuation  and  travel  time  files. 

3.  tlen(j),  j  =  1,  10  -  Time  window  length  used  for  the  spectral  estimate  for  each 
wave. 

4.  gvell(j)  and  gvel2(j),  j  =  1,  10  -  Range  of  group  velocities  used  for  the  spectral 
estimate  for  each  wave  (if  tlen(j)  is  less  than  zero). 

5.  numpthmed  -  Number  of  different  path  media. 

6.  pthmed(i),  i  =  1,  numpthmed  -  Path  media  names  (e.g.,  stable  and  tectonic). 

7.  tbl.atn(i),  i  -  1,  numpthmed  -  File  names  for  attenuation  data  tables. 

8.  tbl.tim(i),  i  =  1,  numpthmed  -  File  names  for  travel  time  data  tables. 

9.  atnstd(j,i),  j  =  1,  10  and  i  =  1,  numpthmed  -  Standard  deviation  of  log  attenua¬ 
tion. 

10.  sorcor(j.i),  j  =  1,  10  and  i  =  1,  numpthmed  -  Log  source  correction  factor  for 
teleseismic  distances. 

11.  reccor(j,i),  j  =  1,  10  and  i  =  1,  numpthmed  -  Log  receiver  correction  factor  for 
teleseismic  distances. 

12.  timsd(j,i),  j  =  L  10  and  i  -  1,  numpthmed  -  Travel  time  standard  deviation  (s). 

13.  azsd(j,i),  j  =  1,  10  and  i  =  1,  numpthmed  -  Azimuth  standard  deviation  (degrees). 

Path  Grid  File  (.pgr) 

1.  title  -  Title  of  the  path  grid  file. 

2.  numpgrlat  numpgrlon  -  Number  of  latitudes  and  longitudes  in  the  path  medium 
grid. 

3.  pgrlatl  pgrionl  -  Southernmost  latitude  and  westernmost  longitude  in  the  path 
medium  grid. 

4.  pgrlatdei  pgrlondel  -  Latitude  increment  and  longitude  increment  for  the  path 
medium  grid. 

5.  pthmedname(l),  l  -  1,  numpgrlat  x  numpgrlon  -  Path  medium  name  (e.g.,  stable 
or  tectonic)  for  each  grid  element. 

6.  pgrlai(i)  and  pgrlon(j),  i  -  1,  numpgrlat  and  j  -  1,  numpgrlon  -  Latitude  and 
longitude  of  the  southwest  comer  of  each  grid  element. 

Attenuation  File  (.atn) 

1.  title  -  Title  of  the  attenuation  file. 

2.  rumatnfrq  -  Number  of  frequencies  in  the  attenuation  file. 

3.  atnfrq(i),  i  -  1,  numatnfrq  -  Frequencies  (Hz). 
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4.  numatndst  -  Number  of  distances  in  the  attenuation  file. 

5.  atndst(j),  j  =  1,  numatndst  -  Distances  (degrees). 

6.  atnamp(j,i),  j  ~  1,  numatndst  and  i  =  1,  numatnfrq  -  Attenuation  amplitude  (log 
m-1). 


Travel  Time  File  (.tim) 

1.  title  -  Tide  of  the  travel  time  file. 

2.  numtimdph  -  Number  of  depths  in  the  travel  time  file. 

3.  timdph(i),  i  =  1,  numtimdph  -  Depths  (km). 

4.  numtimdst  -  Number  of  distances  in  the  travel  time  file. 

5.  timdst(j),  j  =  1,  numtimdst  -  Distances  (degrees). 

6.  tratim(j,i),  j  =  1,  numtimdst  and  i  -  1,  numtimdph  -  Travel  times  (s). 

Site/Station  Data.  The  site/station  data  define  the  station  parameters  (e.g.,  station 
locations,  station  corrections,  station  reliability,  signal-to-noise  ratios  required  for 
detection,  etc.)  and  the  local  site  response  and  array  gain  functions.  There  are  two 
input  file  types  in  this  category:  the  station  file  and  die  site  response  file.  The  station 
file  is  indexed  in  the  control  file.  The  site  response  file  (for  each  station  and  wave)  is 
indexed  in  die  station  file.  The  station  file  also  indexes  the  ambient  noise  and  noise 
factor  files  (next  section).  Tables  A.10-A.11  show  the  format  of  the  site/stauon  input 
files.  The  following  list  defines  each  input  parameter. 

Station  File  (.sta) 

1.  title  -  Tide  of  the  station  file. 

2.  dir  -  Directory  for  the  site  response,  noise,  and  noise  factor  files. 

3.  tbl.nfc(i),  i  -  1, 10  -  File  name  for  the  noise  factor  (coda  decay  rate)  file. 

4.  prewav(i),  i  -  1,  10  ~  Phase  identification  of  the  previous  arrival  (e.g.,  i  =  5  for 
Lg  and  prewav(5)  -  Sn ). 

5.  numsta  -  Number  of  stations  in  the  network. 

6.  staid(j),j  =  1,  numsta  -  Station  identification  (e.g.,  ANTO). 

7.  tbl.noi(j),  j  -  1,  numsta  -  File  name  for  the  ambient  noise  power  spectral  density. 

8.  noisd(j),  j  =  1,  numsta  ~  Standard  deviation  of  the  log  ambient  noise  amplitude. 

9.  stalat(j)  stalon(j),  j  =  /.  numsta  -  Station  latitude  and  longitude. 

10.  starel(j),  j  =  1,  numsta  -  Station  reliability  (<  1.0). 

11.  starho(j)>j  =  1,  numsta  -  Density  of  the  receiver  medium  (kg/m3). 

12.  stacor(ij),  i  -  1, 10  and  j  ~  I,  numsta  -  Station  correction  to  log  amplitude. 

13.  stsigsd(ij),  i  =  1,  10  and  j  =  1,  numsta  -  Station-specific  standard  deviation  of 
the  log  signal  amplitude. 
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Table  A.10.  Station  File 


STATION  FILE 


title 

dir 

tbl.nfc(l)  prewav(l) 


tbl.nfc(lO)  prewav(lO) 
numsta 

staid(  1 )  tbl.noi(  1 )  noisd(  1 ) 
stalat(l)  stalon(l)  starel(l) 
starho(l) 

stacor(l,l)  stsigsd(l.l)  snrth(lj) 

statimsd(l,l)  staazsd(l,i)  stavel(l,l)  tbl.srs(l,l) 


stacor(10,l)  stsigsd(10,l)  snrth(10,l) 

statimsd(10,l)  staazsd(10,l)  stavel(10,l)  tbl.srs(10,l) 


staid( numsta)  tbl.noi( numsta)  noisd( numsta) 
stalat( numsta)  stalonf numsta)  starel( numsta) 
starho(numsta) 

stacor(l .numsta)  stsigsd(  1 , numsta)  snrth(l, numsta) 

statimsdfl, numsta)  staazsd(l, numsta)  stavel(l , numsta)  tbl.srs(l , numsta) 


stacor(  10, numsta)  stsigsdf  10, numsta)  snrthf  10, numsta) 

statimsd(  10, numsta)  staazsd(  10, numsta)  stavel(  10, numsta)  tbl.srs(10,  numsta) 


Table  A.11.  Site  Response  File 


SITE  RESPONSE  FILE 

Suffix:  .srs 


title 

numsrs 

srsfrqf  1 )  srsamp(  1 ) 


srsfrq(  numsrs)  srsamp( numsrs) 
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14.  snrth(ij),  i  =  1,  10  and  j  =  1,  numsta  -  Signal-to-noise  ratio  threshold  for  detec¬ 
tion. 

15.  statimsd(ij),  i  =  1,  10  and  j  =  1,  numsta  -  Station-specific  travel  time  standard 
deviation  (s). 

16  staazsd(ij),  i  =  1,  10  and  j  =  1,  numsta  -  Station-specific  azimuth  standard  devi¬ 
ation  (degrees) 

17.  stavel(ij),  i  =  1,  10  and  j  =  1,  numsta  -  Velocity  of  the  receiver  medium  (in 
m/s). 

18.  tbl.srs(ij),  i  =  1,  10  and  j  =  1,  numsta  -  File  name  for  the  site  response  data  file. 

Site  Response  File  (.srs) 

1.  title  -  Title  of  the  site  response  file. 

2.  numsrs  -  Number  of  frequencies  in  the  site  response  file. 

3.  srsfrq(i),  i  =  1,  numsrs  -  Frequency  (Hz). 

4.  srsamp(i),  i  =  1,  numsrs  -  Log  amplitude  of  the  site  response. 

Noise  Data.  The  noise  data  define  the  ambient  station  noise  spectra  and  coda  decay 
rates  that  determine  the  noise  levels  for  secondary  phases  (see  Section  2).  There  are 
two  input  file  types  in  this  category:  the  noise  file  and  the  noise  factor  file  (coda  decay 
rate).  Both  of  these  files  are  indexed  in  the  station  file.  There  is  a  separate  ambient 
noise  file  for  each  station,  and  a  separate  noise  factor  file  for  each  wave  type.  Tables 
A.12-A.13  show  the  format  of  the  noise  input  files.  The  following  list  defines  each 
input  parameter. 


Ambient  Noise  File  (.noi) 

1.  title  -  Title  of  the  ambient  noise  file. 

2.  numnoi  -  Number  of  frequencies  in  the  ambient  noise  file. 

3.  noifrq(i),  i  =  1,  numnoi  -  Frequency  (Hz) 

4.  noiamp(i),  i  =  1,  numnoi  -  Log  noise  power  spectral  density  (nm2/Hz.) 

Noise  Factor  File  for  Secondary  Phases  (.nfc) 

1.  title  -  Title  of  the  noise  factor  file. 

2.  numnfc  -  Number  of  distances  in  the  noise  factor  file. 

3.  nfodst(i),  i  =  1,  numnfc  -  Distance  (degrees) 

4.  nfqfac(i),  i  =  1,  numnfc  -  Noise  factor  (or  coda  decay  rate). 
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Table  A.12.  Ambient  Noise  File 


AMBIENT  NOISE  FILE 

Suffix:  .noi 

title 

numnoi 

noifrq(l)  noiamp(l) 

noifrq( numnoi)  noiamp( numnoi) 
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Table  A.13.  Noise  Factor  File 


NOISE  FACTOR  FILE 

Suffix:  .nfc 

title 

numnfc 

nfcdst(l)  nfcfac(l) 

nfcdst( numnfc)  nfcfacf numnfc) 
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OUTPUT  DATA 


There  are  two  NetSim  output  files:  a  graphics  file  (suffix  =  ".gra")  and  a  list  file  (suffix 
=  ".1st").  These  suffixes  are  appended  to  the  output  prefix  specified  in  the  control  file. 
The  graphics  file  is  an  ASCII  file  that  is  read  by  SnapCon  to  create  contour  maps  of 
the  NetSim  output  (see  Section  3).  The  list  file  is  also  an  ASCII  file  and  it  includes 
individual  station  probabilities,  the  frequency  of  the  maximum  signal- to-noise  ratio, 
and  location  importances  for  each  wave,  station  and  epicenter. 

Graphics  File.  The  graphics  output  file  contains  information  from  the  control  file  to 
identify  the  run  parameters.  More  importantly,  it  contains  epicenters,  detection  thres¬ 
holds,  probabilities,  and  location  uncertainties  that  are  used  as  input  to  SnapCon  for 
plotting  contour  maps.  Tables  A.  14  shows  the  format  of  the  graphics  output  file. 
Most  of  the  parameters  in  the  ".gra"  file  were  described  above  (control  input).  Only 
parameters  not  described  previously  are  identified  in  the  following  list. 

Graphics  File  (.gra) 

1.  episize(i),  i  =  1,  numepi  -  Event  size  in  units  specified  by  sizetype.  These  are 
detection  thresholds  for  a  NetSim  threshold  run. 

2.  fnetprob(i),  i  -  1,  numepi  -  Probability  of  detection  for  the  network. 

3.  fracloc(i),  i  =  1,  numepi  -  Fraction  of  Monte  Carlo  realizations  with  sufficient 
data  to  constrain  location. 

4.  smajl(i)  and  smaj2(i),  i  -  1,  numepi  -  Length  of  the  semi-major  axis  (km)  for 
prctlo  and  prcthi ,  respectively  (smajl(i)  =  smaj2(i)  for  probability-weighted  loca¬ 
tion  uncertainties). 

5.  sminl(i)  and  smin2(i),  i  =  1,  numepi  -  Length  of  the  semi-minor  axis  (km)  for 
prctlo  and  prcthi ,  respectively  (sminl(i)  =  smin2(i)  for  probability-weighted  loca¬ 
tion  uncertainties). 

6.  strkl(i)  and  strk2(i),  i  -  1,  numepi  -  Strike  of  the  semi-rdnor  axis  for  prctlo  and 
prcthi,  respectively  (strkl(i)  =  strk2(i)  for  probability- weighted  location  uncertain¬ 
ties). 

7.  depthl(i)  and  deptl2(i),  i  =  1,  numepi  -  Depth  uncertainty  (km)  for  prctlo  and 
prcthi,  respectively  (depthl(i)  =  depth2(i )  for  probability-weighted  location  uncer¬ 
tainties). 

List  File.  The  list  output  file  contains  detailed  information  about  the  detection  and 
location  parameters  at  each  station  and  for  each  epicenter.  Tables  A.  15  shows  the  for¬ 
mat  of  the  list  output  file.  The  first  page  lists  the  control  data  for  the  run.  The  second 
page  lists  the  wave-independent  station  data,  and  this  is  followed  by  one  page  of  sta¬ 
tion  data  for  each  wave  used  in  the  NetSim  run.  After  that,  there  are  two  pages  for 
each  epicenter.  The  first  of  these  Usts  for  each  wave  at  the  individual  stations  the  pro¬ 
bability  of  detection,  the  frequency  of  the  maximum  signal-to-noise  ratio,  and  the 
amplitude  standard  deviation.  The  second  page  for  each  epicenter  lists  for  each  wave 
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Table  A.14.  Graphics  Hie 


Title  of  Run 

Date  and  Time  of  Run 

Control  Hie  Name 

Epicenter  Hie  Name 

Source  Media  File  Name 

Source  Grid  File  Name 

Path  Media  File  Name 

Path  Grid  File  Name 

Station  File  Name 

Output  Prefix 

Output  Directory 

Input  Directory 

Run  Type 

Sub-run  Type 

Detection  Criteria 

Waves  to  Use  in  Location 

Event  Size  Type 

Event  Size 

Minimum  Event  Size 

Maximum  Event  Size 

Probability  of  Detection 

Low  Probability  Cutoff 

Regional  Distance 

Ijocaiion  Confidence  Level 

Number  of  Monte  Carlo  Simulations 

Low  %  Monte  Carlo 

High  %  Monte  Carlo 

Fixed  or  Free  Depth 

eventdepth 

Number  of  Frequencies 
Run  Frequencies 


numsta 

slaid(l)  stalat(l)  slalon(J) 


GRAPHICS  FILE 
Suffix:  .gra 

:  t  itle 

:  idate,  itime 
:  filtjcnt 
:  file  j pi 
:  filejmd 
■•filtjgr 
:  filtpmd 
-.filtpgr 
:  filejta 
:  outpre 
:  outdir 
:  dir 

:  runlype 
:  subtype 
:  detcrit 
■  loewaves 
:  sizetype 
:  evsise 
:  evmin 
:  evmax 
:  probdet 
:  eutlo 
:  regdist 
:  probloc 
:  numcarlo 
:  prctlo 
:  prethi 
:  depthfix 
:  eventdepth 
:  numrmfrq 
■■  runfrq(l) 


runfrq(numrunfrq) 


slaidfnumsla)  stalat(numsta)  stalon(numsta) 
numepi 

epilat(l)  epilon(l)  tpnize(I)  fneiprob{l)  fracloc(l) 
smajl(l)  sminl(l)  strkJ(l)  depthl(l) 
smafi(l)  smin2(l)  Mk2(l)  depth2(l) 


epilat(numepi)  epilon(mimepi)  episize(numepi)  fiutprobfnumepi)  fracloc(numepi) 
smajl(numepi)  sminl(numepi)  strkl (numepi)  depth! (numepi) 
smaj2(numepi)  smin2(numepi)  strk2(numepi)  deplh2(nutnepi) 


Table  A.  15.  List  File 


LIST  FILE  (Page  1) 

.  Suffix •  Jrt 

Title  of  Run 

title 

Dale  and  Time  of  Run 

ielate,  itime 

Control  Hie  Name 

file.au 

Epicenter  File  Name 

file.epi 

Source  Media  File  Name 

filesmd 

Source  Grid  Hie  Name 

fUejgr 

Path  Media  Hie  Name 

file.pmd 

Path  Grid  File  Name 

Station  Hie  Name 

fiesta 

Output  Prefix 

outpre 

Output  Directory 

outdir 

Input  Directory 

dir 

Run  Type 

runtype 

Sub-run  Type 

subtype 

Detection  Criteria 

deterit 

Waves  to  Use  in  Location 

loewaves 

Event  Size  Type 

sizetype 

Event  Size 

evsize 

Minimum  Event  Size 

evmin 

Maximum  Event  Size 

evntax 

Probability  of  Detection 

probdet 

Low  Probability  Cutoff 

cutlo 

Regional  Distance 

regdist 

Location  Confidence  Level 

probloc 

Number  of  Monte  Carlo  Simulations 

numcarlo 

Low  %  Monte  Carlo 

prcllo 

High  %  Monte  Carlo 

prelhi 

Fixed  or  Free  Depth 

depthftx 

evenldepth 

evenldepth 

Number  of  Frequencies 

numrunfrq 

Run  Frequencies 

runfrq(l) 

runfrq(nutnrunfrq) 

LIST  FILE  (Page  2) 

Suffix:  1st 


STATION  DATA 

STATION  FILE  DIRECTORY:  stadir 


I 

IDENT 

LAT. 

LONG. 

RET  TAR 

DENSITY 

NOISE 

SD 

AMBIENT-NOISE 

TABLE 

isia 

staid 

stolon 

starel 

starho 

1  oisd 

tbljtoi 
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LIST  FILE  (Page*  3-X) 
Suffix:  Jst 


STATION  TABLE  (wavename  -WAVE) 


I 

IDENT 

STATION 

CORRECTION 

SI N 
RATIO 

WAVE 

VELOCITY 

SIGNAL 

SD 

TIME 

SD 

AZIMUTH 

SD 

SITE  RESPONSE 
TABLE 

ista 

staid 

stacor 

snrth 

stavel 

stsigsd 

statimsd 

staatsd 

tblsrs 

LIST  FILE  (Page  X+l) 

Suffix:  Jst 


EPICENTER  iepi  LAT=  epilal  LONG=  epilon 

EVENT  SIZE 
NUM.  OF  ITERATIONS 
SEMI-MAJOR  AXIS  (LO) 
SEMI-MINOR  AXIS  (LO) 
STRIKE  (LO) 
DEPTH-INTERVAL  (LO) 
MONTE  CARLO  FRACTION 


ANG. 

STATION 

STATION 

DIST 

AZM 

PROB. 

ista  staid 

stadist 

staaxm 

staprob 

episize 

NETWORK  PROBABILITY  (%) 

:  fnetprob 

number 

ITER.  BOUNDS  EXCEEDED 

:  outbound 

smajto 

SEMI-MAJOR  AXIS  (HI) 

:  smajhi 

sminlo 

SEMI-MINOR  AXIS  (HI) 

:  sminhi 

strklo 

STRIKE  (HI) 

:  strkhi 

depthlo 

fracloc 

DEPTH-INTERVAL  (HI) 

:  depthhi 

wavename  ■ 

WAVE 

FREQ.  OF 

AMPL 

MAX.  S/N 

ATTENSD 

frqmax 

ampsd 

LIST  FILE  (Page  X+2) 
Suffix:  Jst 


EPICENTER  iepi 


LAT=  epilal  LONG=  epilon 


wavename  -WAVE 

ANG. 

EPICENTER 

DEPTH 

STATION 

DIST 

AZM 

A'.'D 

AZSD  IMPORT. 

IMPORT. 

and  tpiimp 


ista  staid 


stadist 


staaxm 


timed 


depthimp 


at  each  station  the  arrival  time  and  azimuth  standard  deviations,  and  the  epicenter  and 
depth  importances  (see  Ciervo  et  al.,  1985). 
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DISTRIBUTION  LIST 
FOR  UNCLASSIFIED  REPORTS 
DARPA-FUNDED  PROJECTS 
(Last  Revised:  26  Sep  89) 


RECIPIENT 


NUMBER  OF  COPIES 


DEPARTMENT  OF  DEFENSE 


DARPA/NMRO  2 

ATTN:  Dr.  R.  Alewine  and  Dr.  R.  Blandford 
1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 
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Washington,  D.C.  20340-6158 
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Shock  Physics  Directorate/SD 
Washington,  D.C.  20305-1000 

Defense  Technical  Information  Center  2 
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Alexandria,  VA  22314 


DEPARTMENT  OF  THE  AIR  FORCE 
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Bldg  410,  Room  C222 

Bolling  AFB,  Washington,  D.C.  20332-6448 

AFTACySIINFO  1 

Patrick  AFB,  FL  32925-6001 

AFTAC/IT  3 

Patrick  AFB,  FL  32925-6001 
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Kirkland  AFB,  NM  87171-6008 
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Washington,  D.C.  20505 


U.S.  Arms  Control  and  Disarmament  Agency 
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Massachusetts  Institute  of  Technology 
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